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ABSTRACT 


This report describes simplified mathematical models of the Space 
Shuttle's Orbital Maneuvering System (OMS), Reaction Control System (UCS), 
and on-orbit Digital Autopilot (DAP) that have been incorporated in the High- 
Fidelity Relative Motion Program (HFRMP) for the HP-9825A desk-top calculator. 
Comparisons are made between data generated by the HFRMP and by the Space 
Shuttle Functional Simulator (SSFS), which models the cited Shuttle systems 
in much greater detail. These data include propellant requirements for 
representative translational maneuvers, rotational maneuvers, and attitude 
maintenance options. Also included are data relating to on-orbit trajectory 
deviations induced by RCS translational cross coupling. Potential close- 
range stationkeeping problems that are suggested by HFRMP simulations of 80 
millisecond (as opposed to 40 millisecond) DAP cycle effects are described. 

The principal function of the HFRMP is to serve as a flight design tool in 
the area of proximity operations. However, the simplified models it uses may 
be of value in other applications, for instance in ground-based orbit deter- 
mination (navigational) calculations. 
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1 . INTRODUCTION 


The HP-9825A High Fidelity Relative Motion Program (HFRMP) is a 12- 
degrees-of-freedom trajectory integrator (6 degrees of freedom for each of 
two vehicles) which generates digital and graphical data to describe the rel- 
ative motion of the Space Shuttle Orbiter and a free-flying payload. These 
data are obtained by differencing the geocentric states of the individual 
vehicles, computed to a numerical precision of 12 decimal digits, with respect 
to an oblate earth whose gravitational model includes the second harmonic co- 
efficient (Jg). The state of the two-vehicle system is computed as a function 
of time by means of a fourth-order Runge-Kutta numerical integration scheme 
which uses quaternions to define both the rotational and the translational 
states of each vehicle (Reference 1). 

The payload is modeled geometrically as a cylinder whose length and diam- 
eter are specified by program input. The Orbiter and the payload are treated 
as rigid bodies whose individual mass properties (gross weight, moments and 

products of inertia, and center of gravity location) are specified by program 

+• 

input and are assumed to remain constant during the HFRMP run. Gravity gra- 
dient torque is included in the rotational equations of motion for both vehicles. 
At the user's option, aerodynamic torque and drag can also be included for 
either or both vehicles. Aerodynamic forces normal to the relative wind vec- 
tor are ignored. 

The atmosphere is assumed to rotate with the earth, and is modeled as 
a function of geodetic altitude by a curve fit of the 1962 Standard density 
profile (Reference 2). The density profile can be modified by a program input 
factor to account for the major effects of solar activity. Aerodynamic drag 
and moment coefficients for the Orbiter are computed as curve-fit functions of 
its attitude with respect to the relative wind vector (References 3-5). 


'An exception is made when integrating the trajectory of an upper stage during 
a major translational maneuver such as an IUS solid rocket motor burn. In 
such a case, the payload gross weight is decremented during the burn to re- 
flect the consumption of propellant. 
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Aerodynamically, the payload is modeled as a flat plate whose size and 
shape are determined by projecting its cylindrical outline onto a plane normal 
to the relative wind vector. The payload drag coefficient is assumed to be 
2.0, based on its projected frontal area. Aerodynamic effects on the Orbiter 
and the payload can be modified (or cancelled entirely) by means of input fac- 
tors which are applied uniformly to all aerodynamic forces and torques that are 
computed for the specified vehicle. 

Several options are provided for defining the initial state of the two- 
vehicle system. The translational state of the Orbiter can be described either 
in terms of osculating orbit elements referenced to the Mean of 1950.0 (M50) geo 
centric equatorial frame, or in terms of invariant orbit elements (Reference 6) 
measured in the Mean of Launch Date (MLD) equatorial frame. The initial atti- 
tude of the Orbiter is defined by pitch, yaw, and roll angles (taken in that 
order) referenced to the rotating Shuttle-centered local -vertical (SLV) coor- 
dinate system. The Orbiter's angular velocity,, measured in terms of rate 
components about its body axes, can be defined relative to either the M50 
(inertial) frame or the SLV (rotating) frame. 

The initial translational state of the payload is defined by rectangular 
position and velocity components which are measured relative to the Orbiter's 
center of gravity (CG). At the user's option, these components can be measured 
in the SLV coordinate system, or in the Shuttle body (SBY) coordinate system. 

The payload's initial pitch, y aw, and roll angles can be referenced either to 
the payload-centered local-vertical ( PL V ) system, or to the SBY system. The 
payload's angular rate components about its body axes can be defined relative 
to the M50, the PLV, or the SBY frame. 

The initial state of the system can be advanced through up to 40 flight 
profile segments, each of arbitrary length, which are defined by the user at 

4 * 

the beginning of the HFRMP run. At the beginning of any segment the user may 


+ A1 1 input data, including the flight profile definition, are saved in disk 
files whence they can be recalled (and edited, if necessary ) for use in 
subsequent runs. 
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command the application of an impulsive (i.e,, instantaneous) increment to the 

.t. 

angular rate of either or both venicles. 1 In this regard, the user may specify 
a particular rate increment (INCR), a desired rate with respect to inertial 
space (IR), or a desired rate with respect to the local -vertical frame of the 
vehicle in question (LVR). In all cases, the components of the desired rate 
or rate increment are measured about the body axes of the vehicle in question. 

After applying the specified angular velocity increment (if any) at the 
beginning of the flight segment, the HFRMP then (for each vehicle independently 
of the other) 

1. allows the attitude to drift (D) under the influence of inertia and 
natural torques, or 

2. performs inertial rate hold (IRH) control (i.e., maintains a constant 
angular velocity relative to inertial space), or 

3. performs local -vertical rate hold (LVRH) control (i.e., maintains a 
constant angular velocity relative to the rotating 1 ocal -verti cal 
frame of the vehicle being controlled) 

for the duration of the segment, depending on the user's specifications for 
that segment. 

When the IRH or the LVRH attitude-maintenance option is specified for the 
Orbiter, a simplified RCS/DAP model is used to compute average values for the 
propellant consumption rates and translational cross-coupling accelerations that 
result from the intermittent thruster firings which are required to apply the 
necessary control torques. The model takes into account the mass properties 
of the Orbiter, the electrical width (an integer multiple of the DAP cycle time) 
and the effective width (the duration of steady-state acceleration) of the RCS 
thruster pulses, and the width of the attitude deadband about each of the Orbiter's 


’’’Commanded angular velocity impulses, and the linear velocity impulses which 
they may induce as a result of RCS translational cross coupling, are the only 
types of state variable discontinuity that are permitted by the HFRMP. These 
are allowed only at the beginning of a flight profile segment. 
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body axes. Deadbands can be changed from segment to segment in the flight 
profile, as can the selection of primary or vernier thrusters and the mode 
of cross-coupling compensation. Translational cross-coupling accelerations 
are integrated along with those produced by gravity, aerodynamics , and com- 
manded translational thrust. They are reflected in the output data by the 
flight path deviations they produce. Propellant consumption rates are also 
integrated (but not subtracted from the Orbiter gross weight), and the accum- 
ulated expenditures are tabulated, along with other data, a user-specified 
time intervals. Propellant consumption is broken down according to source 
(forward, aft left, or aft right tank) and function (translational or rota- 
tional control). 

When the IRH or the LVRH attitude-maintenance option is specified for the 
payload, the magnitudes of the necessary control torques are integrated and the 
accumulated rotational impulse (measured in pound-foot-seconds) in the positive 
and negative direction about each body axis is printed along with the other 
digital output data. Since no specific method of implementation is modeled, it 
is not possible to compute propellant consumption rates or cross-coupling effects 
that may result from payload attitude control . 

Translational thrust acceleration of either or both vehicles can be 
commanded at the beginning of any f 1 i ght profile segment. Payload translational 
thrust is always applied in the direction of the payload’s +X body axis and 
is assumed to be directed through the CG. Once initiated, payload thrust ac- 
celeration continues until all of the rocket motor propellant is consumed, as 
determined by a table of flow rates versus burn time. 

Translational acceleration of the Orbiter is initiated by commanding 
ignition of either or both of the OMS engines (L, R, or L+R) and/or by firing 
primary RCS thrusters to produce thrust 'nominally in the positive or negative 
directions of the Orbiter body axes (+X, -X, +Y, -Y, +2, or -Z). Once initiated, 
Orbiter translational acceleration is applied continuously at the nominal steady- 
state level, throughout the duration of the flight profile segment. Detailed 
descriptions of the OMS and RCS models are contained in Sections 2 and 3. 
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2. OMS MODEL 


The thrust vectors of the OMS engines are assumed always to pass through 
their respective gimbal centers (STA 1518, BL +88, WL 492) as defined by Ref- 
erence 7 in terms of Orbiter station coordinates. When a single-engine OMS 
burn is commanded, the specified engine is gimbal ed to direct its thrust vector 
through the Orbiter CG. When both engines are commanded to fire, they are 
gimbaled such that their two thrust vectors will be parallel and have moment 
arms about the CG of equal magnitude and opposite direction. Thus, in all cases, 
the OMS engines produce zero net torque. 

Throughout the total duration of any flight profile segment in which it 
is commanded to fire, each OMS engine is assumed to deliver its nominal steady- 
state thrust magnitude of 6000 pounds, and to consume propellant at its nominal 
rate of 19.16 pounds/second (Reference 7). Flow rates are integrated and the 
accumulated propellant expenditure is listed separately for each engine in the 
HFRMP output data. However, the propellant expenditure is not subtracted from 
the Orbiter's gross weight. Therefore, this simplified model will not generate 
accurate data for a very large OMS burn where the propellant expenditure repre- 
sents a significant fraction of the Orbiter gross weight. 
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3. RCS MODEL 


Table 1 contains the basic data which are used to compute forces and 
torques produced by individual RCS jets. Each row in the table contains data 
for a particular jet. The first column on the left contains index numbers 
which are used internally by the HFRMP. The second column contains the jet 
identification mnemonics which are defined in Reference 7. The next three 
columns contain steady-state thrust components parallel to the Orbiter body 
axes. These are followed by three columns which contain the station coordi- 
nates of a thrust application point (i.e., a point on the jet's line of action). 
These coordinates are used, in conjunction with the quantity in the last col- 
umn, to calculate torque about the Orbiter CG. 

The data in Table 1 were derived from those shown in Tables 2 and 3, 
which are reproduced from References 7 and 8. The thrust components in Table 1 
represent a summation of the basic thrust components from Table 2 and the in- 
cremental thrust components (due to plume impingement on Orbiter surfaces) 
from Table 3. For each jet listed in Table 3, the thrust application point 
coordinates in Table 1 were obtained by the following procedure: 

1. Using the force components and application points from Table 2, a 
basic force moment vector was calculated about the reference CG 
location (STA 1076.7, BL 0, WL 374.1) which is associated with 
Table 3. 

2. The torque increment due to impingement (from Table 3) was added to 
the torque calculated in step 1. 

3. The resultant torque vector was resolved into two components; one 
parallel to the total force vector (columns 3-5 of Table 1), and 
the other normal to the total force vector. 

4. The parallel torque component was divided by the magnitude of the total 
force vector. The result, which has the dimensions of feet, appears 

in the rightmost column of Table 1. 

5. A line of action (effective thrust line) was calculated for the total 

force vector such that its moment about the reference CG was equal to 
the normal component of the total torque vector. The station coordi- 
nates of the new thrust application point (columns b-8 of Table 1 ) 
represent the point on the calculated thrust line that is closest to 
the thrust application point originally designated in Table 2. ^ 
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Table 1. RCS Thruster Data 


THRUSTER 

NO. 

THRUSTER 

ID 

mm 

an 

WKBM 

STA 

(IN) 

BL 

(IN) 

WL 

(IN) 

c 

(FT) 

1 

F2F 

-879.4 

-26.2 

119.9 

306.72 

14 . 

65 

392.96 

O.OGOO 

2 

F3F 

-879. 5 

0.0 

122. 7 

306.72 

0. 

00 

394.45 

0.0000 

3 

F1F 

-879.4 

26.2 

119. 9 

306.72 

-14. 

65 

392. 96 

o.oooo 

4 

F1L 

-26.3 

873.6 

18.2 

362.67 

-69. 

50 

373.73 

o.ooou 

5 

F3L 

-21.0 

870.3 

0.5 

364.71 

-71. 

65 

i 359.25 

0 .OOUO 

6 

F2R 

-26.3 

-873.6 

18.2 

362.67 

69. 

50 

1 373.73 

o.ouou 

7 

F4R 

-21.0 

-870.3 

0.5 

364.71 

71. 

65 

i 359.25 

U. UU uu 

a 

F2u 

-32.3 

-11.7 

874.4 

35U.93 

14 . 

39 

413.46 

0.0U00 

9 

F3Q 

-31.9 

0.0 

873.5 

350.92 

0. 

00 

414.53 

U . UOuu 

10 

FlU 

-32.3 

11.7 

874.4 

350.93 

-14. 

39 

413.4b 

0. UUU0 

11 

F2U 


-616.4 

-639.5 

333.84 

61. 

42 

356. 95 

0 . uUUu 

12 

FID 


616.4 

-639. 5 

333.84 

-61. 

42 

35b . y5 

U.ouuu 

13 

F4D 

-24.8 

-612 .6 

-639.4 

348.44; 

66 . 

23 

358.44 

O.OuOu 

14 

F3D 

-24.8 

612.6 

-639.4 

348.44 

-66. 

23 

358.44 

O.OuOU 

15 

R3A 

856.8 

0.0 

151.1 

155 5.29: 

137. 

00 

473.06 

o.ooou 

16 

RlA 

856.8 

0.0 

151.1 

1555.29: 

124 . 

00 

473.06 

O.OuUo 

17 

L3A 

856.8 

0.0 

151.1 

1555.29 

-137. 

00 

473.06 

o.uuoo 

IB 

LlA 

856.8 

0.0 

151.1 

1555. 29 

-124 . 

00 

473.06 

O.OUUU 

19 

L4L 

0.0 

870.5 

-8.4 

1516.06 

-149. 

83 

455.21 

-0.5887 

20 

121 

0.0 

870.5 

-8.4 

1529.07 

-149. 

83 

455.21 

-0.6061 

21 

131 

0.0 

870.5 

-8.4 

1542.07 

-149. 

83 

455.21 

-0.6235 

22 

L1L 

0.0 

870.5 

-8.4 

1555.07 

-149. 

83 

455.21 

-0.6410 

23 

R4R 

0.0 

-870.5 

-8.4 

1516.06 

j 149. 

83 

455.21 

0. 5887 

24 

R2R 

0.0 

-870.5 

-8.4 

1529.07 

! 149. 

83 

455.21 

0.6061 

25 

R3R 

0.0 

-870.5 

-8.4 

1542.07 

; 149. 

83 

455.21 

0.6235 

26 

RlR 

0.0 

-870.5 

-8.4 

1555. 07 

149. 

83 

455.21 

0. 641U 

27 

. 14U 

29.0 

72.0 

870.0 

1520.04 

I -116. 

51 

481.65 

-0.4615 

28 

12U 

29.0 

72.0 

870.0 

1532.96 

-116. 

54 

481.65 

-0.3725 

29 

LIU 

29.0 

72.0 

870.0 

1545187 

-116. 

58 

481.65 

- 0 . 28 36 

30 

R4U 

29.0 

-72.0 

870.0; 

1520.04 

116. 

51 

481.65 

0 .4615 

31 

R2U 

29.0 

-72.0 

870.0, 

1532.96 

116. 

54 

481.65 

0 . 3725 

32 

RlU 

29.0 

-72.0 

870.0 

i 1545.87 

116. 

58 

481.65 

0.283b 

33 

14D 

312.4 

346.8 

-545.7 

14 98.11 

. -101. 

47 

42U.49 

1.7413 

34 

12D 

312.4 

34 6.8 

-543. 7 

1513.68 

-100. 

61 

WBEm 

1.4807 

35 

1 3D 

312.4 

346.8 

-545.7 

1529.23 

-99 

79 

BrrirH 

1.220b 

36 

K4D 

312.4 

-346.8 

-545.7 

14 98.11 

101. 

47 

MBB 

-1. 7413 

37 

R2D 

312.4 

-346.6 

-545.7 

1513.68 

100. 

61 

424.63 

-1.48U7 

38 

R3D 

312.4 

-346.8 

-545.7 

1529.23 

99. 

79 

428 . 76 

-1. 22U8 

39 

F5R 

-0.8 

-17.0 

-17.6 

324.35 

59. 

70 

350.12 

O.UUOO 

40 

F51 

-0.8 

17.0 

-17.6 

324.35 

-59. 

7 U 

350 . 12 

0 . Uu 00 

41 

R5K 

0.U 

-24.0 

-0.6 

1565.00 

14 9. 

8 7 

459.00 

O.OuUO 

42 

151 

0.0 

24.0 

-0. 6 

1565. 0U 

-149. 

87 

459. 00 

, 0 .0 0 U U 

43 

R5D 

0.0 

0.0 

-24.0 

1565.00 

118 

00 

455.44 

O.OUUU 

44 

15D 

0.0 

0.0 

-24. 0 

15b5 . 0 0 

-118. 

0O 

455.44 

0. UUUU 











Table 2. Basic RCS Thrust Data 

(Without Plume Impingement) 


THRUSTER 

ID 

THRUST COMPONENTS, LB 

RESULTANT 

STATION COORDINATES 
THRUST APPLICATION P 

OF 

OINT (IN) 

Fx 

F y 


THRUST, LB 

STA 

BL 

WL 

F2F 

-879.4 

-26.2 

119.9 

887.9 

306.72 

14.65 

392.96 

F3F 

-879.5 

0.0 

122.7 

888.0 

306.72 

0.0 

394.45 

FI F 

-B79.4 

26.2 

119.9 

B87.9 

306.72 

-14.65 

392.96 

F1L 

-26.3 

873.6 

18.2 

874.2 

362.67 

-69.50 

373.73 

F3L 

-21.0 

870.3 

0.5 

870.6 

364.71 

i -71.65 

359.25 

F2R 

-26.3 

-873.6 

18.2 

874.2 

362.67 

69.50 

373.73 

F4R 

-21.0 

-870.3 

0.5 

870.6 

364.71 

71.65 

359.25 

F2U 

-32.3 

-11.7 

874.4 

875.1 

350.93 

14.39 

413.46 

F3U 

-31.9 

0.0 

873.5 

874.1 

350.92 

0.0 

414.53 

F1U 

-32.3 

11.7 

874.4 

875.1 

350.93 

-14.39 

413.46 

F2D 

-28.0 

-616.4 

-639.5 

888.6 

333.84 

61.42 

356.95 

FID 

-28.0 

616.4 

-639.5 

888.6 

333.84 

-61.42 

356.95 

FAD 

-24.8 

-612.6 

-639.4 

885.9 

348.44 

66.23 

358.44 

F3D 

-24.8 

612.6 

-639.4 

885.9 

348.44 

-66.23 

358.44 

F5F. 

-0.8 

-17.0 

-17.6 

24.5 

324.35 

59.70 

350.12 

F5L 

-0.8 

17.0 

-17.6 

24.5 

324.35 

-59.70 

350.12 

R3A 

856. 8 


151.1 

870.0 

1555.29 

137.00 

473.06 

RJ A 

856.8 

0.0 

151.1 

870.0 

1555.29 

124.00 

473.06 

L3A 

B56 . 8 

0.0 

151.1 

B70.0 

1555.29 

-137.00 

473.06 

L1A 

856.8 

0.0 

151.1 

870.0 

15551.29 

-124.00 

473.06 

L4L 

0.0 

870.5 

-22.4 

870.8 

1516.00 

-149. B7 

459.00 

L2L 

0.0 

870.5 

-22.4 

870.8 

1529 .00 

-149.87 

459.00 

L3L 

0.0 

870.5 

-22.4 

870.8 

1542j.00 

-149.87 

459.00 

L1L 

0.0 

870.5 

-22.4 

870.8 

1555>00 

-149.87 

459 .00 

R4K 

0.0 

-870.5 

-22.4 

870. E 

1516.00 

149.87 

459.00 

R2R 

0.0 

-870.5 

-22.4 

870.8 

1529.00 

149.87 

459.00 

R3R 

0.0 

-670.5 

-22.4 

870.6 

1542.00 

149.87 

459.00 

R1R 

0.0 

-870.5 

-22.4 

870.8 

1555.00 

149.87 

459.00 

L4U 

0.0 

0.0 

670.0 

870.0 

1516.00 

-132.00 

4B0.50 

L2U 

0.0 

0.0 

870.0 

870.0 

1529 L 00 

-132.00 

480.50 

LID 

0.0 

0.0 

870.0 

870.0 

1542.00 

-132.00 

480.50 


0.0 

0.0 

870.0 

870.0 

1516 >00 

132.00 

480.50 

R2U 

0.0 

0.0 

870.0 

870.0 

1529.00 

132. 00 

480.50 

RIO 

0.0 

0.0 

870.0 

870.0 

1542.00 

132.00 

480.50 

LAD 

170.4 

291. B 

-801.7 

870.0 

1516.00 

-111.95 

437.40 

L2D 

170.4 

291.8 

-801.7 

870.0 

1529.00 

- 111.00 

440.00 

L3D 

170.4 

291.8 

-801.7 

870.0 

1542.00 

-110.06 

442.60 

R4D 

170.4 

-291.8 

-801.7 

870.0 

1526.00 

111.95 

437.40 

! R2D 

170.4 

-291.8 

-801.7 

870.0 

1529.00 

111.00 

440 . 00 

| R3D 

170.4 

-291.8 

-801.7 

870.0 

1542.00 

110.06 

442.60 

! L5D 

0.0 

0.0 

-24.0 

24.0 

1565,00 

-118.00 

455.44 

R5D 

0.0 

0.0 

-24.0 

24.0 

1565.00 

118.00 

455.44 

L5L 

0.0 

24 .0 

-0.6 

24.0 

1565.00 

-149.67 

459.00 

F.5R 

0.0 

-24.0 

-0 .6 

. 24.0 

1565.00 

149.87 

459.00 


(Reproduced from Reference 7) 
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Table 3. Force and Moment Increments Due to Plume Impingement 
(CG (3 STA 1076.7, BL 0, WL 374.1 ) 


THRUSTER 

i ID 

1 

A 

FORCE (LB) 

A MOMENT (LB-FT) 


F 

y 

F z 


S' 

L 

z 

1 



14 

-450 

iSgljg 





14 

-450 



L3L 



14 

-450 

■ 1 


LI L 



14 

-450 



L4U 

29 

72 


1755 

II 

-2780 

L2U 

29 

72 


1755 

mmm 

-2780 

LIU 

29 

72 


1755 

IMi& 

-2700 

L4D 

142 

55 

256 

-2520 


-1395 

L2D 

142 

55 

256 

-2520 

wmm 

-1395 

L3D 

142 

55 

256 

-2520 

10,480 

-1395 

R4R 



14 

450 



R2R 



14 

450 



R3R 



14 

450 



R1R 



14 

450 



R4U 

29 

-72 


-1755 


2780 

R2U 

29 

-72 


-1755 


2780 

R1U 

29 

-72 


-1755 


2780 

R4D 

142 

-55 

256 

2.520 

1 0 ,480 

1395 

R2D 

142 

-55 

256 

2520 

10,480 

1395 

R3D 

142 

-55 

256 

2520 

10,480 

1395 


(Reproduced from Reference 8) 
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For a CG location defined by arbitrary station coordinates (STA c g, BL C g, 

WL ) , the torque produced by a particular jet j is computed from the equation 
*-9 


l j ‘ r j x f j ♦ Vj • 


where 


V F yj 


•(STA. - STA )/l 2 
v 3 eg 

(BL. - BL )/ 1 2 
v 3 eg" 

•(WL. - WL )/12 
v J eg" 


and where the values of F . , F ., F . , STA., BL., WL . , and C- are obtained from 

XJ y J Z J J J J j 

columns 3-9 of Table 1. 
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4. DAP MODEL 


4.1 JET SELECT TABLES 

For each flight profile segment, the HFRMP user must specify which one 
of three basic combinations of thrusters is to be used for attitude and/or 
RCS translational control of the Orbiter. The available options are designated 
V (vernier jets), P (primary jets), and PZI (primary jets with +Z thrusters 
inhibited). Corresponding to each of those options is a jet-select table 
(Tables 4-6) which identifies the particular jet or combination of jets that is 
to be fired in response to each of the six possible translational acceleration 
commands (+X, -X, +Y, -Y, +Z, -Z) and the six possible rotational acceleration 
commands (+R0L, -ROL, +PCH, -PCH , +YAW, -YAW). Jets are identified by the 
mnemonics listed in the second column of Table 1. The jet select tables are 
not routinely available for modification by the user. However, an editing rou- 
tine is available in the HFRMP software system, and changes can be made (or 
additional options can be provided) with little difficulty. 

As indicated in Table 4 by the absence of any jet designations for the 
execution of translational acceleration commands, the V option (vernier jets) 
can be used only for attitude control. The P option (Table 5) and the PZI 
option (Table 6) can be used for translational and/or rotational control. 

In the PZI option, no jets are fired that would expel propellant directly 
upward with respect to the Orbiter body. Translational acceleration in the down 
ward direction, if commanded, is achieved (at a comparatively high propellant 
cost) by firing the +X and -X thrusters simultaneously. The cant angles of 
the tX jet and -X jet thrust lines produce a small net acceleration in the +Z 
(downward) direction. This option normally is used only when the Orbiter is 
maneuvering in the near vicinity of a payload that must be protected from jet 
plume impingement. 

4.2 RESPONSE MATRICES 

During HFRMP trajectory integration, steady-state accelerations and 
propellant consumption rates that result from RCS translational and rotational 
control commands are obtained from a pre-calculated "response matrix". As 
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Table 4. Vernier (V) Jet Select Table 


1 

CMD 

i 



THRUSTERS 

TO 

BE FIRED 



1 

2 

3 

Km 

5 

6 

7 

8 

1 

+X 
-X 
+i 
-£ 
i + z 
: -z 

; +Rujj, 

! -RJiii, 

i fRCri 
i-PCn. 
+ £ AW 
-if AW 

LOU 

k5jj 

r‘3K 

i/5D 

tOK 

L5L 

t'ib 

k5jj 








ORIGINAL PAGE IS 
OF POOR QUALITY, 
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ORIGINAL PACK'S 

Q E EOOR OUAU.V 


Table 5. Primary (P) Jet Select Table 


! CMD 

THRUSTERS TO BE FIRED | 

1 

2 

3 

4 

5 

6 

7 

8 

+x 

RlA 

Lii A 

m m 






1 -x 

f 2 f 

c‘ij’ 







1 + ¥ 

t'lL 

; i-j 4 Li 

‘ tl 






j -Y 

£' 2 R 

r4k 

mam 






1 + 2 

PJu 

Li4u 







1 ~' L 

t'lu 

r' 2 L; 

■Mill 

l 2 D 

R 4 J 

k 2 u 



j +ROL 

li 4 u 

R 4 U 







1 -ROL 

L 4 U 

R 4 u 







1 + PCH 

c'lb 

C'ZU 

Li 4 U 

R 4 U 





-PCH 

PJu 

L 4 u 

R 4 l) 






+ YAW 

PlL 

K 4 R 







-YAtf 

j?' 2 k 

J_j 4 Li 
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Table 6. Primary with +Z Thrusters Inhibited (PZI) 
Jet Select Table 


rMn 

THRUSTERS TO BE FIRED ~] 

UmU - 

1 

2 

3 

4 

5 

6 

7 

8 

[ 

i +x 

KlA 

LilA 







-A 

R2R 

RlR 








RiL 

ii 4-L* 







-Y 

c7h 

K 4 K 







:i 

R2H 

Rlr 

i<liA 

lil A 





-L 

c'iu 

r 1 2D 

L 4U 

L 2u 

K4d 

KZU 



+ROi-i 

j_i 4U 








j-RJii 

k4u 








+*»Cu 

flu 

R2d 







r-RCii 

J_1 4 Li 

rt4D 







+ ZAX 

r’lL 

i<4t< 







-YA rt 

R2k 

Li 4L 
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indicated in Tables 7-14 for a typical set of Orbiter mass properties which are 
listed in Table 15, the response to a particular command is conditioned by the 
jet select table (V, P, or PZI), and it is further conditioned by the user- 
specified mode of cross-coupling compensation. The available compensation modes 
are NONE (no compensation), ROT (rotational compensation only), and FULL (rota- 
tional and translational compensation). 

The FULL compensation mode is not applicable when the vernier jets are 

4 . 

selected; therefore, eight different response matrices are sufficient to sup- 
port all of the Orbiter control options that are currently provided in the HFRMP. 
Only one of these eight matrices is required to support trajectory integration 
in a given flight profile segment. AIT eight matrices are stored in disk files 
prior to the initiation of trajectory integration, and the appropriate one is 
read into calculator memory at the beginning of each segment. All eight matrices 
are re-calculated whenever the user makes any kind of change to the Shuttle Data 
Base (where the Orbiter mass properties are stored). 

The response matrices are 12 x 12 arrays of real numbers. In a given 
matrix, each row corresponds to a particular translational or rotational con- 
trol command. The first three columns of each matrix contain the body-axis 
components of the steady-state linear acceleration vector, and columns 4-6 
contain the corresponding components of angular acceleration. Columns 7-9 
contain rates of flow from the forward, aft left, and aft right tanks (in that 
order), for propellant consumption that is chargeable to translational control. 
Columns 10-12 contain similar rates for propellant consumption that is charge- 
able to rotational control . 

4.2.1 Uncompensated Response Calculations 

Linear accelerations appearing in the uncompensated response matrices 
(Tables 7-9) are calculated by use of the equation 


4 - 

The compensation mode is internally defaulted to ROT if the user specifies 
FULL compensation for the vernier jets. 

4 * 4 * 

In Tables 7-14, each matrix is partitioned into two 12 x 6 arrays, simply 
for the purpose of fitting it conveniently on a page. Command identifica- 
tion mnemonics are tabulated in a column to the left of each matrix parti- 
tion. 
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Table 7. Uncompensated Response Matrix for Jet Select Option V 


CMD 

LINEAR ACCELERATION (i 

r T/SEC 2 ) 

ANGULAR ACCELERATION (] 

RAD/SEC 2 ) 

a x 

a y 

a z 

“x 

a y 

a z 


0 .UU 00 0 

0.00000 

0. OOOOO 

0.00000 

0. uOOOO 

0. uOuou 

-x 

0 .OOOOU 

0 . 0 0000 

0. 00000 

0. 00000 

O.OUOOO 

0. OOoOU 

+Y 

0 . 0 0 0 0 0 

0.00000 

0. OOOOO 

0.00000 

0. 00000 

0. OoOOu 

-I 

U. 00000 

O.UOOOO 

0. OOOOO 

0. uOOOU 

O.OUOOO 

0. OOUuu 

; ■ 

0.00000 

0. 00000 

0.00000 

0. oOuOO 

0. uOuUU 

0 . OOoOU 


U . 00 00 0 

O.uOGGG 

0 . oouoo 

0. UOOOU 

U. UOOOU 

0 . U 0 0 u o 

+ KU1. 

o.uuouO 

O.ouOOO 

— 0.0 0336 

U . OOu^.7 

-0 .00015 

u. uuoOl 

-RUii 

, 0 . 0 0 0 u u 

0. OOOUO 

-0.00306 

— U . 00 J2 1 

-0.00 ui 5 

—0. uOOul 

+ir*Cii 

~U . U002o 

0 .oOOOO 

-u. OUbob 

0 . uuo 00 

0.00 03 5 

u. Uuuuu 

-FCii 

0. 00 000 

0 . 0 0 0 0 0 

-0.00772 

— U .OoOUu 

-0 . 00 U2o 

0.0 OUOO 

+ ¥A»V 

U . 00 00 0 

— 0 . oOido 

-o . ooolo 

- 0 . U u U 1 0 

-0. UuOOO 

0.00 Ul3 

-YA.V 

0 . 0000 0 

o . O03o6 

-0 . uGUiO 

U. oUUlb 

— 0. uOuOu 

— U . 0 U Ol 



PROPELLANT CONSUMPTION RATES (LB/SEC) | 

CMD 

TRANSLATIONAL CONTROL 

ROTATIONAL CONTROL | 


FWD 

AFT LEFT 

AFT RIGHT 

FWD 

AFT LEFT 

AFT RIGHT 


TANK 

TANK 

TANK 

TANK 

TANK 

TANK 

+A. 

0 . OOOOO 

o . u 0 G u U 

U.OuuOu 

u. uoooo 

0. uuoUu 

u. uuo uu 

-x 

0 .00000 

0 . oouoo 

O.OUOOO 

0. OOOOO 

O.UUOUO 

0. OOOOU 

+r 

0 .00000 

0 . uoooo 

0 . uoooo 

0.00000 

0.00000 

0 . 0 u 0 0 0 

-t 

0.00000 

0. OOOOO 

o. uoooo 

0 . ououo 

0. OOOUO 

U. OUUOU 

+ 2 

0.00000 

0.00000 

0 . oouoo 

0.00000 

0.00000 

0. OOUOO 

-2 

0.00000 

0 . uoooo 

0. OOOOU 

0.00000 

0.00000 

o . o UG oo 

+RUJL 

0 .00 00 0 

0.00000 

U. OOOOO 

0.00000 

0.09 230 

u . Oouuu 

-RuL 

0.00000 

0.00000 

0. OOOOO 

0.00000 

0. OOOOU 

0.09230 

+UCri 

0.00000 

0.00000 

0. UUUGO 

0.13460 

0.00000 

0. OOOOO 

-ROii 

0.00000 

0. UOOOO 

0. OOOOO 

0.00000 

0.09230 

0.09230 

+ 3f A M 

U .00000 

0. OOOOU 

0. UOUOU 

U.UOuoO 

0. OOOuO 

0.09230 

-1 A/f 

0 .00000 

0.00000 

0. UOOOO 

0. OOOOO 

0.09230 

o.uouuu 



(See Table 15 for Orbiter Mass Properties) 
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Table 8. Uncompensated Response Matrix for Jet Select Option P 


CMD 

LINEAR ACCELERATION (FT/SEC 2 ) | 

1 ANGULAR ACCELERATION (RAD/SEC 2 ) 

a x 

a y 

a z 

a x 

a y 

a Z 

+ X 

0 . 2 75u4 

0 . U 0 0 0 0 

U . U 4 d o 1 

-u . uooul 

-0 . u 0u J 3 

0 . 0 0 0 ol 

-x 

-0. 28 291 

O.UOOOO 

0. 03357 

-o.ouoo<; 

-0 . UU 2 ll 

-a , oouoi 

+ Y 

-0.00423 

0 . 2 d 0 5 5 

0. 0U158 

o . o o 7 o 7 

-Q.U0030 

0.00365 

-Y 

-J . UU423 

-0 . 2d ud5 

0. 00153 

-U . 00703 

-U.00U31 

-0.00365 

+ 2 

0.00420 

0.00000 

0. 4 2040 

-0. 00007 

U.uulOl 

-0. OUOuU 

-2 

0.19200 

0.00000 

—0.55635 

0. OU010 

Q.OUUol 

0. 00 00 1 

+ROL 

0 . 0 5 4 9 2 

0 . 04 420 

0. 03217 

0.01591 

0 .00184 

-0. 0 0 050 

-RGii 

0 . 05492 

-0.04420 

0.05217 

-0. 01590 

0. 00131 

0. 00 050 

+PCd 

0.00032 

0 . 0 0 0 0 0 

0. 07415 

0 . 0 u 0 1 3 

0.02 225 

0. UOoOu 

-PCH 

0. 09537 

O.UOOOO 

— 0. U3505 

— 0 . OOuOd 

-0. U1437 

U. OUOuU 

+ YAA 

-0.UG423 

0.0 UUdo 

0. 00153 

— 0 . 0 03 45 

-0.00031 

0. 01230 

-YA/i 

-0 . u04 23 

-0. 00050 

0.00153 

u. U U 3 4 5 

-u.uUOju 

. _ _ . 

— 0. 01233 



j PROPELLANT CONSUMPTION RATES (LB/SEC) 

CMD 

TRANSLATIONAL CONTROL 

ROTATIONAL CONTROL 


FWD 

AFT LEFT 

AFT RIGHT 

FWD 

AFT LEFT 

AFT RIGHT 


TANK 

TANK 

TANK 

TANK 

TANK 

TANK 

+ X 

U „ u 0 0 0 u 

3. lu 71u 

3. lu71U 

0 . 0 U 000 

J . 0 0 0 0 0 

U . U OuUU 

-x 

6. 21420 

u . 0 0 0 0 0 

u. UdOuO 

0 . OOUOO 

0. uOuuU 

0 . 0 Oo uO 

+Y 

d. 107x0 

3. 10/lu 

o. uuUuu 

0.00 0 uO 

U. ooOOu 

0 . U 0 U U 0 

-i 

3.10 710 

O.UOOUO 

3. 10 / lo 

0. OoO Ou 

0 . UOuOo 

0 . U Ou uo 


3. 10710 

3. 10 /lu 

3. 10 710 

0. UUUOO 

0 . 0 J 0 0 0 

0. UOuuU 

-2 

0.21x20 

o . 2 14 20 

6 . 2142 u 

0 . 0 U 0 0 0 

0.00000 

0. uuOOu 

+RJi. 

O.OO00O 

0 . 0 0 0 0 U 

0. UOOOO 

J. UUOUO 

j. 1 o 7 1 o 

3.10710 

-ROL 

0 . OOdU 0 

u.uoouu 

0. UUO00 

0. OUUOu 

3. 10 710 

3 . lu71U 

+PCR 

0 . uOOOO 

0 . u 0 0 0 0 

0. uOuou 

6 . 21 42 u 

3. lO 71 0 

3 . lo 710 

-PCil 

0.00000 

u . 0 Ul) 00 

0.00000 

3. 10710 

3 . lu 71 0 

3. 10710 


0. 00000 

0.00000 

0. OOUOO 

3. 10710 

0 . 0 0 0 U 0 

3.10710 

-*AW 

0. 00000 

o. ooooo 

O.UOOOO 

3.10710 

3. 10710 

i 

0 . 0 0 0 o 0 


(See Table 15 for Orbiter Mass Properties) 



























Table 9. Uncompensated Response Matrix for Jet Select Option PZI 


CMD 

Linear acceleration ( 

FT/SEC 2 ) 

ANGULAR ACCELERATION (RAD/SEC 2 ) | 

a x 

a y 

a Z 

a x 

a y 

a z 

+x 

0 . 27 584 

U.uOOOO 

0. 04dbl 

-0 . UOU Ul 

— 0.0U033 

0. UOUul 

-x 

-0 . 2bzy 1 

0 . UOUuO 

0. 03857 

-0. 0 00 02 

-U.UUZll 

— 0 . 0 uUol 

+Y 

- 0 . u 0 4 2 3 

0.28055 

0 . 00 15d 

0. uu707 

- 0 . 0 0 0 3 0 

0. uo3ab 

-Y 

— u . 0 U 4 23 

-u . zd ub5 

0.00158 

-0. 0o7o8 

-U. udool 

-0 . 0o3b5 

+.4 

-0 . Ji) 727 

0 . uOUUO 

0. Ud 71 8 

— u . uuOu3 

-U.U 

- 0 . oouoo 

-6 

u .xyZOO 

0.00000 

— U . do uii 

0. OuUlU 

0. OuUol 

0 . 0 0 UOi 

+AvJXi 

0 . 0r>0 25 

0 . ud57cJ 

- 0 . 0 d 7 7 b 

U . 00 id 6 7 

-0 . 0UZ34 

— u » u u 12 3 

— RuL 

0 . UuU25 

- 0 . 0 05 7 o 

-0 . Ud 778 

- 0 . u 06 8 8 

- 0 . 0 U 2 'J o 

0. Uul23 

+RCti 

-0.00301 

0 . OU000 

-0. 205 74 

0. OUUlz 

0 . U 12 u 3 

U. UuuUU 

-ir'Ctl 

U. lU05u 

U.UOOOd 

-0. 175do 

-0. OUUUl 

-0.00530 

0 . OOUOO 

+ YAW 

-0 . 004 23 

0. 00050 

0. 00158 

-0. 00345 

-0. 00031 

0. 012 3b 

— YAW 

-J .U0423 

-0 . OOOuu 

0.00158 

O.OU345 

-u.uuOlu 

-0. 01238 


CMD 

PROPELLANT CONSUMPTION RATES (LB/SEC) | 

TRANSLATIONAL CONTROL 

ROTATIONAL CONTROL | 

FWD 

TANK 

AFT LEFT 
TANK 

AFT RIGHT 
TANK 

FWD 

TANK 

AFT LEFT 
TANK 

AFT RIGHT 
TANK 

+x 

0 .uOOuO 

3.10710 

3. 10 710 

0. 00000 

0.00000 

. u. Oouoo 

-X 

6.21420 

U. 00000 

0. UOQOQ 

0. 00000 

0. UOOUO 

0. UOUUO 

+Y 

3.10710 

3.10710 

0 . 0 0 0 0 0 

0 . U 0 0 Uu 

0. 00000 

0 . oouoo 

-Y 

3.10710 

0. uOUUO 

j. 1 0 7 1 u 

0. UOUUO 

0. OOODu 

0. UUOOU 

+ Z 

6 . 21420 

3. 10710 

3. 10710 

0. UUOOu 

0.00000 

0. OOOOU 

■67 

a 

6 . 21420 

6 . 214 20 

•o . 21420 

0. UUOOU 

O.UUOOO 

o. uOoou 

+xiuL 

0 .80000 

O.UUOOO 

0. OOUOO 

0. UUOOU 

3. iu 71o 

0. UUUuO 

— Ruii 

0. uGOOO 

0.00000 

U. OOUOO 

u. ooouu 

u. uu UOU 

3 . 1 0 7 1 U 

+-PCii 

0. OOUOO | 

U. 00000 

0. OU00U 

6. 2142U 

0 . 0 0 ou u 

u. oooou 

-KJti 

0 . 0 0 U u 0 

u. OOUOO 

U. uUUUu 

0 . OOuOO 

3. Iu71u 

3.107x0 

+ YAW 

O.uUOOO 

0. UOUOU 

0. UOOUO 

3.10710 

u. UOOUO 

3.1O710 

-YAW 

U. UOOUO 

0.00000 

O.UUUGO 

3. 1U71U 

3.10710 

0 . u UOuO 


(See Table 15 for Orbiter Mass Properties) 
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Table 10. Rotationally Compensated Response Matrix for Jet Select Option V 


CMD 

LINEAR ACCELERATION ( 

r T/SEC 2 ) 

ANGULAR ACCELERATION ( 

RAD/SEC 2 ) 

a x 

a y 

a z 

a x 

“y 

“z 

+X 
-x 
+ Y 
-1 
+ L 
-L 

+HuiJj 

-Kuij 

+-fc'Oii 

— Jr-wti 

-iaW 

0 .00000 
0 .00000 
0 . 00000 
U. 00000 
0.0 0000 
0 .UGOOO 
-U. 00011 
-0. uOUll 
-0 . J 0 U 2 o 
U. UUOU 0 
-U .OOOOo 
— 0 .00 00 6 

0.00000 
0.000U0 
0.00000 
0.00000 
0. O0000 
0.00000 
0 . 00 023 

-0 .00026 

— u.oOouu 
0 . 0 0000 

— U . u 0336 
u . U U 3 8 o 

0. 00000 
0.00000 
0. OOOOO 
0. ooOOO 
0. OOOOO 
0. OOOOO 

— 0 • uO o22 

— O.uUozJ 
-0.00571 
-u. uO 77 3 
-0 . uU 37 2 
-0.00375 

0. OOOOO 
U. OOOOO 
0. OOOOO 
0. 000 JO 
0. OOOOO 
u . oOooo 

u . 00020 

- U . 0 0 0 2 0 

— u . u UO 0 0 
0. uOGOO 

-0 . u UU UU 

0 . oooou 

0. OUOOO 
0. UU uOU 
0. 00 00 0 
U . 00 00 0 
0 . 0 0 U 0 u 
0. OUOOO 
0. Oouo 0 
0 . 00 0 0 0 
0. UO u3 5 
-u . 00 0 2y 
0 . OU uO 0 

0 . uouoo 

u . uu ouu 

0. OOUUO 

0 , U uOUu 
0 . 0 0 U 0 J 
u. UuuOO 
0. uO jOU 
0. OOUuO 
0 . u U u U 0 
U. UUuUU 
0 . 0 U 0 0 0 
0. 0u0l4 
-0. uUUl4 


■ '1 

CMD 

PROPELLANT CONSUMPTION RATES (LB/SEC) j 

TRANSLATIONAL. CONTROL 

ROTATIONAL CONTROL | 

FWD 

TANK j 

AFT LEFT 
TANK 

AFT RIGHT 
TANK 

FWD 

TANK 

AFT LEFT 
TANK 

AFT RIGHT 
TANK 

+ A 

0.00000 

U.UOOOO 

0. UOOOU 

0. UUUUU 

U. UUUUU 

U . UUuUU 

-X 

0 . UO OOU 

U.UOOOO 

0. UOUUU 

u.uuuOu 

0. UUUUU 

U. UUUUU 

+7 

0 .00000 

O.OUOUO 

0. uUOOO 

0. JOOUO 

0 . 00.00 0 

O.OuOOu 

-L 

0 . 0 00 0 0 

0 . oooou 

0. OOOOO 

0. OOUOU 

0 . UOUu u 

U. UOUUU 


0.00000 

0.00000 

0. UUOOO 

0 . OOOOO 

0 . uu uu u 

0. OUOUU 


0 .00000 

0.00000 

0. OOOOO 

0. OOOOO 

0. OOOOO 

0. OOUuO 

+RGL 

0 .0000 U 

0. UOUOO 

U.OOUOO 

0. 07632 

U. U3914 

O.OOOUO 

-RuL. 

0. OOOOO 

0 . OOOOO 

0. OOOOO 

0.07708 

0. OOOOO 

0. U99od 

+ ±»CH 

0.00000 

0.00000 

0. UOOOO 

0. 18460 

-0 .00004 

0.00109 

-Jk>CH 

0.00000 

0.00000 

0. OOOOO 

0. UUOOO 

0. 09240 

0. 09230 

+YAtf 

0.00000 

0.00000 

U.UOOOO 

0.04600 

0,05296 

0,09230 

i -tm 

0.00000 

0. OOOOO 

0 . UUOOO 

0. 04635 

0.05230 

0.05341 


(See Table 15 for Orbiter Mass Properties) 
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Table 11. Rotational ly Compensated Response Matrix for Jet Select Option P 


CMD 

LINEAR ACCELERATION ( 

FT/SEC 2 ) 

ANGULAR ACCELERATION (RAD/SEC 2 ) | 

a x 

a y 

a z 

“x 

a y 

a z 

+x 

0.27570 

0.O0001 

0. 04973 

0 . 00000 

0,. UUOO0 

0. oouuu 

-x 

-0. 20274 

0.00003 

0.04572 

-U.OOO00 

— 0.00000 

0. oooou 

+ Y 

0.02797 

0.25772 

0. 02 775 

0.00000 

0.00000 

0. UUUOO 

-Y 

0.02811 

-0.25769 

0. 02780 

-0.00000 

0.00000 

0. 0 0 U 0 0 

+ Z 

0.01172 

0.00021 

0. 41798 

0.00000 

-u.ouoou 

0. 00000 

-z 

0.19652 

-0.00029 

-0.55801 

0. 00000 

-0.00000 

0.00000 

+ROL 

0.06689 

0.04422 

0. 04777 

0.01576 

0 . 0O000 

-0. 00OUO 

-ROL 

0.06670 

-0.04422 

0. 04783 

-0.01578 

-0.00000 

0. 00000 

+PCH 

0.00077 

-0.00036 

0. 07458 

0.00000 

0.02226 

— 0. OOOuO 

-PCti 

0.09567 

0.00024 

-0.03477 

0.00000 

-0. 01436 

-0. 00000 

+YAW 

0.00768 

0. 01009 

0.01258 

-0.00000 

-0.00000 

0.01227 

-YAW 

0 .00767 

- 0 . 0 100 8 

0. 01256 

0.00000 

0. GO OU 0 

-0. 01227 


ro 

o 



PROPELLANT CONSUMPTION RATES (LB/SEC) | 

CMD 

TRANSLATIONAL CONTROL 

ROTATIONAL CONTROL f 


FWD 

AFT LEFT 

AFT RIGHT 

FWD 

AFT LEFT 

AFT RIGHT 


TANK 

TANK 

TANK 

TANK 

TANK 

TANK 

+x 

0.00000 

3.10710 

3. 10710 

0. 09260 

0. 04962 

0.04812 

-X 

6.21420 

0.00000 

0 . O 0000 

0. 59174 

0.30238 

0. 304 22 

+Y 

3.10710 

3.10710 

0. 00000 

1. 22951 

2 . 7 a oUl 

1.80583 

-Y 

3.10710 

0. 00000 

3.10710 

1. 23580 

1. 8l3o6 

2 . 79308 

+z 

3. 107.10 

3. 10710 

3. 10 710 

0 . 23769 

0 . 251zl 

0 . 25242 

— z 

6.21420 

o. 21420 

b. 21420 

0 . 18739 

0 . Id 548 

U. 15659 

+ROL 

0.00000 

0.00000 

0. 00000 

0 . 52082 

3. 5 OZ'l 4 

3.0X752 

-ROl 

0 . 00 00 U 

0.00000 

0. 00000 

0. 51593 

3.o 22 5 3 

3 . 4 9 o 7 2 

+PCH 

0.00000 

0.00000 

0. 00000 

6. 21561 

3. 13390 

3 . 13249 

-PCri 

0.00 00 0 

0.00000 

0 . UOOOU 

3. 10673 

3. 12 33 z 

3.12369 

PYAW 

0.00000 

O.oOOOO 

0 . O 000 U 

3. 08110 

0 . 06 II 8 

3 . 76823 

-YAW 

O.OOOOO 

0.00000 

O.OOOOO 

3. 08020 

3. 76692 

0. b5982 


(See Table 15 for Orbiter Mass Properties) 
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Table 12. Rotationally Compensated Response Matrix for Jet Select Option PZI 



1 TNFAR AFTFI FRATT ON ( 

FT/SEC 2 ) 

i AMGIII A R Arm FRATT ON / RAD/FFf.^l 1 

CMD 

ll_l 11 » < 




■ » * 


a x 

a y 

a z 

a X 

a y 

a z 

+X 

U .27551 

0 . 00003 

0. 04306 

-0 . 000 00 

— 0. OOOOO 

-o. ooOuo 

-x 

-U . 2o4l0 

U. 0U0U4 

0. U034o 

-0 .00000 

—0. uuoou 

0. OoOuO 

+ £ 

U . U o 3 d i 

0.2100a 

-0. 17340 

0. 000 00 

0. 00 oou 

0. OoOoO 

-i 

0 . Ud320 

-0 . 2114 7 

-0.17220 

-u. 00000 

-0. ououo 

0. Ooooo 

+z 

-0 . OUaol 

0 .00007 

U. 0 4 65b< 

-0 . 0U00O 

0. UOOOU 

o. o uooo 

- L 

U . 2U320 

- U . u 0 05 1 

-0.57725 

0 . 00000 

-U. OuOuj 

0. OOUOO 

+RuL 

0 .04772 

0 . 0 35 d 3 

-0 . I35a2 

1 0 . U 0 035 

-0. OOOOO 

0. OOuOU 

-ROli 

U. 04 771 

— 0 . 0 5 5d 3 

- 0 . 13 o U 1 

—0 . OUobl 

-0 . OUU 00 

0 . OOuoO 

+ RCii 

-J . uo all 

-0 . uolol 

-0.20731 

0 . OOOOO 

0. 012 b 4 

o . 0 o u o u 

-FCii 

0. 10054 

0.00004 

-0 . 175ui 

0. UudOO 

-0.00550 

0 . U U 0 0 0 

+ *A,J 

U . UlObJ 

0 . U2 335 

-0. 07113 

- 0 . oouoo 

0 . OOuOO 

0 . Ull7o 

■ —£Aii 

o.ulyai 

— J . U 2 o b u 1 

— 0. 0 7157 

0.0 0 0 0 U 

0 . 0 0 UO u 

-0. oll7a 


PROPELLANT CONSUMPTION RATES (LB/SEC) 


CMD 

TRANSLATIONAL CONTROL 

ROTATIONAL CONTROL 


FWD 

AFT LEFT 

AFT RIGHT 

FWD 

AFT LEFT 

AFT RIGHT 


TANK 

TANK 

TANK 

TANK 

TANK 

TANK 

+x 

0 . 0 0 o 0 0 

3. 10 710 ! 

3. 10 /lu 

u . ±0443 

0 • U \J D u 

0 . U U 2 o / 

-X 

6 . 21420 

0 . 0 0 0 0 0 

0. UOOOU 

1. 04 6 5 7 

0. 01244 

0. U 12.2o 

+ £ 

3 • 1 u 7 1 0 

3 .1 0 7 1 u 

- u. ooooo 

3 . 3 o 7 z y 

1.30124 

3 . a b 3 v a 

-£ 

3 . 10710 

0 . 0 0 0 0 0 

3 . 1 0 / 1 0 

3 . 2 J 0 a 1 

3. 33737 

1. 2 0330 

+' L 

o . 21420 

3 . 1 0 7 1 0 

3. 10710 

1.20723 

u . 01530 

0. Ol25u 

-2 

d . 214 20 

o . 2 14 2 u 

•j . 214 2 0 

0.03107 

U . 5 2 0 u 5 

o . 3 7 2 a 1 

+RJL 

0 . UO JO J 

0.00000 

U. OOOOO 

1 . 7 o 4 o 5 

3.10710 

0 . 3 u3a o 

- dxj Li 

0 . uu 00 0 

J . 0 0 u 0 0 

u . OooOO 

1. 77134 

0 . 3 J j J J 

3.10710 

+±”30 

0 . 0 0 0 0 0 

U. uOOOO 

0. OOuOO 

6 . 2 2 u a 3 

0 . 0 U 0 1 3 

0. U 3 5 5 5 

-RC3 

0.0 00 00 

u . 0 0 0 0 0 

U. UOOOO 

U. uOuol 

3 . luo57 

3 . 1 o 7 1 u 

+ £ Aw 

0 . 0 0 0 0 0 

o . o u o o o 

0 . u J i) 0 0 

3.5 7734 

1. 5a344 

3. 10 /lu 

—£h W 

0. uOOOO 

o . o oooo 

0.0 000 0 

3 . y 3 1 1 a 

3.10710 

1.30457 


(See Table 15 for Orbi ter Mass Properties) 

























Table 13, Fully Compensated Response Matrix for Jet Select Option P 


ro 

ro 


P 0 

^ g 

S§ 

ss 

o ,1; 

i- >*“ 

P a 

H 

t< ss 


CMD 

LINEAR ACCELERATION ( 

FT/SEC 2 ) 

ANGULAR ACCELERATION (RAD/SEC 2 ) 

a x 

a y 

a Z 

“x 

a y 

a z 

+ X 

0 . 2:9 3 3 1 

U . U 0 0 0 0 

o.oouoo 

0. OOUGO 

-0. OuuoO 

o.oouoo 

-X 

— 0 . 2 3 :6 5 4 

-0 .00000 

-u . ooooo 

o. ooooo 

: - o . u o o o o 

0. 00 00 0 

+y 

0 . uOOO 0 

0.25770 

— 0.00 0 00 

0 . 0 0 0 0 0 

0 . ooooo 

O.OCUOu 

—t 

U. 00000 

-0.257 70 

-0.00000 

-0. 0 00 00 

o.oooou 

0. 000 JJ 


0 . 0 u 0 0 0 

o. uuoou 

0.41590 

0 . 0 0 J 0 0 

0 . 00 00 0 

0. UO JuU 


-0 . 00000 

0.00000 

-0 . 5 26 20 

o.oouoo 

-O.OOUOO 

U . 0 0 u 0 0 

+KJL. 

-0.00000 

0.O00U0 

—0 . UUOOU 

0. 01576 

0 . 00 0 0 0 

- u . o 0 0 0 y 

-A -JL 

-0 . 0 0 JOG 

- u . 0 0 0 0 0 

-0. UOOJO 

— 0 . Ol 5 7 d 

-u . 00 J 0 0 

0 . UOOud 

+RCu 

-0 . OOUO 0 

0. 00000 

-0.00000 

0 . Ou UOU 

U. 0222b 

- 0 . 0 J U U 0 

~L J Cn 

0. JOOOO 

0.00000 

— 0 .00 Go 0 

0 . OOUUU 

— 0 . ul4 Jo 

— 0. UOuOu 

+ 2A4 

0 .uOOOO 

0 . u 00 00 

0 . 00 000 

— O.uOOUO 

0 . 00 0 J 0 

0. u 12 2 7 

— Zn.V 

-0.00000 

0.O0000 

U.OUOOU 

u. uO JOu 

0 . 00 UU u 

— 0 . 012 2 7 


CMD 

PROPELLANT CONSUMPTION RATES (LB/SEC) j 

TRANSLATIONAL CONTROL 

ROTATIONAL CONTROL 

FWD 

TANK 

AFT LEFT 
TANK 

AFT RIGHT 
TANK 

FWD 

TANK 

AFT LEFT 
TANK 

AFT RIGHT 
TANK 

+A 

0 . 5 5 3 7 u 

3 . o 6 u 9 J 

3 . o b 'J d 0 

u . Iu4d 8 

0 . 06375 

0. 06193 

-X 

6 . 7 2 3 4 6 

0 . 5oyl4 

u . o o y 2 o 

0 . a 080 9 

0 .3155 2 

0.31716 

+ X 

4 . 3 6 7 o 8 

3. 48822 

0. J8ll2 

1.22173 

2 . 838 oO 

1. 8 584 9 

-1 

4.37231 

u . 3819 7 

3. 4890? 

1. J2d34 

1 . 86584 

2 . 84 3 95 

+ L 

3.3 6 7 o 7 

8.10710 

3. 109b4 

0 . 2632 7 

U. 26524 

0. 2o 733 

-2 

10.53718 

6 . 21738 

6. 21420 

0 . 55050 

0.37253 

0. 3o9 92 

+RO jL 

3.39343 

0.75765 

1 . 2 9 J 5 9 

u . y5004 

3 . 935 8y 

4 . 22944 

— RJL 

3 . 3 8y 4 1 

1.29llo 

0.75770 

0. 9431b 

4 .22245 

3. 92362 

+RCxi 

1.51900 

0 . 83 758 

0.38273 

6. 2972 5 

3. lyisy 

3. 18344 

-FCri 

2.6 9 513 

0 . 44032 

0.44342 

3 . 3 4o 12 

3 . 26 0 20 

3. 26223 

+/A7 

0 .61690 

0,17831 

0. 2yy38 

3 . 1 6 3 6 5 

0.75217 

3 . 8 9 7 6 3 

-im 

0 .61564 

0 . 299 43 

0.17787 

3. 1623d 

3 . 3 '9 5 3 9 

0. 75045 


(See Table 15 for Orbiter Mass Properties} 




Table 14. Fully Compensated Response Matrix for Jet Select Option PZI 



- ; . 


-T/SEC 2 ) 

ANGULAR ACCELERATION ( 

RAD/SEC 2 ) 

CMD 

LINEAR ACCELE RATI ON (T 

a x 

a y 

a z 

a x 

“y 

a z 

+ X 

0.29067 

-u. 00000 

0. 00000 

-0. uouou 

-0. UOUOO 

0. uOUuO 

-X 

-0.26288 

0.00000 

u.oouoo 

-0. 00000 

-0.00000 

0. uuuuo 

+Y 

0 .00000 

0.2x120 

-0.00000 

0. QOuOO 

U. 00000 

0.00000 

-Y 

-0.00000 

-0,21120 

-u.oouoo 

-0.00000 

-0.00000 

0 . oouoo 

+ Z 

-0.00000 

-0.00000 

0.04787 

-0.00000 

-0.00000 

- 0 . oouoo 

~'L 

0.00000 

-0.00000 

-0.57553 

0. 00000 

-0.00000 

0.00000 

+ROL 

-0.00000 

-0.00000 

o. oouoo 

0.00655 

-0.00000 

- 0 . ouooo 

-RUL: 

0.00000 

-0.00000 

0. 00000 

-0.00651 

-0.00000 

0.00000 

+PCH 

-0.00000 

-0.00000 

0.00000 

-0.00000 

0.01264 

0.00000 

-PCti 

-0.00000 

-0.00000 

0. 00000 

-0.00000 

-0.00590 

0.00000 

+YAW 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

-0.00000 

0.01176 

-YAW ; 

-0.00000 

-0.00000 

0.00000 

0.00000 

0.00000 

-0.01175 


: CMD 


PROPELLANT CONSUMPTION RATES (LB/SEC) 

1 

TRANSLATIONAL CONTROL 

ROTATIONAL CONTROL | 

FWD 

TANK 

AFT LEFT 
TANK 

AFT RIGHT 
TANK 

FWD 

TANK 

AFT LEFT 
TANK 

AFT RIGHT i 
TANK 

+x 

0.46330 

3.57088 

3.57040 

0.16626 

0.02999 

0.03030 

-x 

6.25160 

0.03693 

0.03740 

1. 04765 

0. 01498 

0.01467 

+ Y 

26.67336 

14.65366 

11.54656 

7.87348 

1. 36162 

3. 91333 

-Y 

26.50964 

11. 46686 

14.57396 

7.82500 

3.89733 

1.34563 


6.21328 

3.20395 

3.20503 

1 . 213 55 

0.01751 

0.01304 

-z 

10 . 6 7678 

6.22715 

6.21420 

0. 79452 

0.34402 

0. 39780 

+RGL 

25. 69667 

12.12785 

12.95239 

7.45309 

4.18962 

0. 70408 

-RUL' 

25.76410 

12.98288 

12.16461 

7.47731 

0.72213 

4.18116 

+ PCri 

27.05628 

14.04470 

14.03455 

11.51797 

0.08930 

0.13194 

-PCti 

25.02923 

11. 76977 

11.77451 

4. 82657 

3. 18031 

3. 16041 

+ YAW 

13.2319 9 

6. 30440 

6. 72368 

6.90517 

2.10436 

3. 308 Jb 

-YAW 

13.34197 

6.77696 

6.35791 

6.93498 

3.31869 

2. 11534 


(See Table 15 for Orbiter Mass Properties) 





Table 15. Typical Orbiter Mass Properties 


ITEM 

DESCRIPTION 

VALUE 

i 

cl 1 villi’ .................................. 

2uuol7.oo0u id 

2 

JL AX • ••••••••• •••«•••••••■• 9 •••••• • • •••• • 

tSci/302.U0 00 Uiiu a— c‘‘i' 2 

3 

Li‘l „ 

p56o37/.U000 obUvj- r'T 2 


lii.. .. .. . . . ..... .... .. .. ... 

O 0 9 4 3 o 7 . u 0 U 0 ljLuj- r‘ T 2 

'j 

ii il 

- 9 71 .0000 SLuvj— ET z 
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a = (32.1 74/ W ) 


where W is the Orbiter gross weight and where eF • , E F . , and EF • represent 

x j yj z j 

the sums of thrust components (from Table 1) for all of the thrusters that are 
designated in the appropriate jet select table (Tables 4-6) to be fired for the 
purpose of implementing the command under consideration. The corresponding 
angular accelerations are given by 

(5) 


( 6 ) 


(7) 


and where eL. represents the sum of the individual torque vectors produced 
J 

by the designated thrusters,, as calculated from Equations (1) through (3). 

The symbols I ,1 , and I represent the Orbiter's principal moments of 

xx yy z z 

inertia, and primed vectors represent those whose components are measured 
along the correspond! - ng principal axes. 

Expressions (5) and (7) are coordinate transformation equations (see 
Reference 1) in which q~ represents the unit quaternion that defines the angular 
displacement of the Orbiter's principal axes of inertia with respect to to 
body axes, and q represents its conjugate. The quaternion q is computed 
internally by the HFRMP in such a manner that the associated rotation of 
coordinates will transform the inertia tensor 


a ~ q o a o q 


where 


V’xx 


L Ay 

L A Z J 


L H 


f J 


= q 0 EL. 0 q 
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SF xj 

sF yj 
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(which is composed of user-specified moments and products of inertia about the 
Orbiter's body axes) into the diagonal form 

[U = 

Propellant flow rates (columns 7-12) are assumed to be 3.1071 Ib/sec for 
each active primary jet and 0.0923 Ib/sec for each active vernier jet. These 
rates are based on the nominal vacuum thrust magnitudes (870 lb and 24 lb) 
and specific impulses (280 sec and 260 sec) that are given in Reference 7. 

Each thruster is assumed always to be fed from its nominal source (tank), no 
provisions are made for simulating propellant cross-feed. 

4.2.2 Compensated Response Calculations 

The compensated response matrices (Tables 10-14) are obtained by per- 
forming elementary row operations on the uncompensated response matrices 
(Tables 7-9). Specifically, appropriate rows in the uncompensated matrix are 
multiplied by positive factors and added to other rows in such a manner as to 
null the unwanted cross-coupling accelerations. The positive factors represent 
idealized ratios (calculated without regard to minimum-impulse or other limi- 
tations) of the firing time required of the compensating jets with respect 
to the firing time of the jets activated by the primary command. 

The results of the computations described in the preceding paragraph 
are such that propellant consumed in the process of nulling cross-coupled 
angular accelerations is charged to rotational control, and that which is 
consumed in the process of nulling cross-coupled linear accel erations is charged 
to translational control, regardless of the nature of the primary command. 


xx 


yy 
o i 


zz 


(9) 



4.3 ATTITUDE CONTROL LOGIC 


For the purpose of (1) applying an angular velocity impulse at the be- 
ginning of a flight profile segment or (2) exercising rate-hold control (IRH 
or LVRH) within a segment, the HFRMP always applies rotational cross-coupling 
compensation, whether or not such compensation has been specified by the user. 1 " 

For the aforementioned purposes, the compensation mode is internally defaulted 
to ROT if the user has specified NONE; otherwise, the user's specification is 
appl ied. 

4.3.1 Angular Velocity Impulses 

At the beginning of each flight profile segment, the HFRMP determines 
the direction and magnitude of the instantaneous angular rate increment that is 
to be applied about each of the Orbiter's body axes. These three increments 
may be supplied directly by the user (the most common specification being zero 
for all three components), or they may be calculated internally to yield user- 
specified final (post-impulse) rates with respect to inertial space or the 
Orbiter's local vertical frame. For each axis, an effective thruster firing 
time is then calculated by dividing an angular acceleration into the angular 
velocity increment, the acceleration being obtained from the appropriate row 
(+R0L, -ROL, +PCH, -PCH, +YAW, or -YAW) and column (4, 5, or 6) of the applicable 
response matrix (Table 10, 11, 12,13, or 14). The entire row of the matrix 
is then multiplied by the effective firing time, and the resulting linear velocity, 
angular velocity, and propellant consumption increments are added immediately 
to the appropriate state variables and propellant accounts. 

4.3.2 Rate-Hold Control 

The HFRMP trajectory integration step size typically is several orders 
of magnitude greater than the DAP cycle time; therefore, it is not practicable 
to simulate the individual thruster firing pulses that are commanded by the 


T The reasoning here is that it makes little difference in the long run whether 
the effects of rotational cross-coupling are removed immediately upon imple- 
mentation of a primary command, or a short time later when the pilot or the 
DAP determines (on the basis of angular velocity and/or displacement errors 
detected by the TMU) that a corrective command is required. 
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real-world DAP for the purpose of maintaining the Orbiter attitude within 
specified deadbands. The long-term effects of pulsed thruster firings are 
simulated in the HFRMP by integrating the average accelerations and propel- 
lant consumption rates that result from such pulses. An idealized and very 
greatly simplified control logic is used for this purpose. The rationale is 
that if the idealized logic is designed to minimize the average propellant 
consumption, subject to the basic constraints of the real system, it will pro- 
duce essentially the same long-term effects as the real-world DAP. 

At a given instant of time, suppose that inertial effects and natural 
torque (due to aerodynamics and the gravity gradient) are such that, if not 
counteracted by the attitude control jets, they would produce an undesired 
angular acceleration component of magnitude v about one of the Orbiter's 
body axes. When the user specifies the IRH option, v represents the appro- 
priate body-axis component of the Orbiter's total natural acceleration (rela- 
tive to inertial space); when LVRH is specified, it represents the difference 
between the natural acceleration of the Orbiter and that of its 1 ocal -vertical 
reference frame. 

For the purpose of calculating average effects of the control system 
response, we treat v as if it were a constant, and write the equations 

U) = tOg + vt .(10) 

and 

0 = 0g + Wgt + , (11 ) 

to describe the uncontrolled angular velocity and angular displacement (about 
the axis under consideration) as functions of time. 

The purpose of the attitude control logic is to maintain the angular 
displacement 0 within the limits of -D and +D, where D represents the user- 
specified deadband width, in such a manner as to minimize the RCS propellant 
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4* 

consumption. Figure i illustrates an idealized case in which 0 q = D and 
“0 * -^“opf where 

-opt = ^ < 12 > 

represents the magnitude of a uniform corrective rate increment that is applied 
at a regular interval 


T max 


4/577 


(13) 


Equations (12) and (13) are obtained by solving Equations (10) and (11) 
for the angular velocity increment that minimizes the average propellant con- 
sumption rate; which (for any value of -2w 0 = Aw <. 4v/Dv" ) is given by 


W=(T/x)W a , (14) 

where W is the steady-state consumption rate of the thrusters which oppose 
the natural acceleration (i.e., is the sum of columns 7-12 from the appropriate 
row of Table 10, IT, 12, 13, or 14), T is the electrical pulse width of the 
corrective thruster firing, and 

t = Aw/ v (15) 

is the interval between pulses. The magnitude of the angular velocity incre- 
ment is related to T by the equation 

Aw = a 7* , (16) 

where 


'The idealization consists of using angular velocity impulses to represent 
thruster firings, and of making the magnitude of wo precisely half the magnitude 
of the corrective impulse Aw. This is the average magnitude of wq, which is 
randomly distributed between 0 and Aw in the real world. 
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T* = T - .014 seconds 


07 ) 


is the effective pulse width and where a is the magnitude of the steady- 
state angular acceleration produced by the control jets that oppose the 
natural acceleration. 

When (15), (16), and (17) are substituted into Equation (14), the 
average propellant consumption rate is found to be related to Aw by the ex- 
pression 

• o*1 n W V 

W = (1 + ) -2- . (13) 

Equation (18) shows clearly that the average propellant consumption rate in- 
creases when the angular displacement is undercontrolled (i.e., when Aw < 

4\/lDv ), as illustrated in Figure 2. If -2wg = Aw > 4 v / r Dv~, then it becomes 
necessary to apply control pulses at both the upper and the lower deadband limits, 
as illustrated in Figure 3. This represents an overcontrolled situation that 
is even more inefficient than undercontrol , because it results in the alternate 
firing of thrusters that oppose each other. Clearly then, in the idealized 
case under consideration, the optimum value of Aw is given by Equation (12), from 
which we obtain 

T Q pt = 4 \f Dv /a + .014 seconds (19) 

for the optimum electrical pulse width of the average corrective impulse. 

So far, we have treated T as if it were a continuous variable. How- 
ever, in the real world, I must be an integer multiple of the DAP cycle 
time C. Therefore, we write 

T = NC , (20) 

where 

N “MAX (1, INT (T opt /C)) : y:- ^ (21 
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Figure 2. Undercontrolled Angular Displacement History 



and where T Q t is obtained from Equation (19). The function MAX (X,Y) repre- 
sents the maximum value of the arguments X and Y, and the function INT (Z) 
represents the integer part of the argument Z. 

Equations (20) and (21) imply that some degree of undercontrol or over- 
control is unavoidably required in the general case. The degree of overcon- 
trol (if any) is reflected in the value of £ , which represents the average 
acceleration magnitude produced by the reinforcing jets (those which produce 
torque of the same sign as the natural acceleration), and which is given by 
one or the other of the two equations 

[If INT (T opt /C) = 0]: £ = (ctT*)3/16D - v (22a) 

and 

[ If INT (T op t /c ) £ °J : I. " 0 , (22b) 

where T* is given by Equation (17). Equation (22a) can be obtained by sub- 
stituting aT* for Aw Q p t and (v + £) for v in Equation (12). The value of £ 
represents the idealized' minimum amount by which the natural acceleration 
v must be supplemented in order to balance the counteractive control pulses 
of fixed magnitude aT* - a (NC - 0.14), 

After having calculated £ , the average magnitude of the counteractive 
jet acceleration is obtained from the equation 

£ = v> + £ , (23) 

which describes the aforementioned balance of average accelerations. 


‘To arrive at Equation (22a) , the acceleration from the reinforcing jets is 
treated as if it acted continuously at the average level, A greater value, 
approaching £ = (aT*) 2 /OQ - v as v approaches zero, would be obtained if £ 
were treated""as the average resulting from uniform angular velocity impulses 
applied at the lower deadband limit. The use of Equation (22a) is justified 
on the basis of the better agreement between the results it produces and those 
produced, by more detailed simulations of the actual onboard DAP logic (Section 
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To calculate the average linear and angular acceleration vector com- 
ponents that result from controlling angular motion about the body axis under 
consideration, the data in columns 1-6 of the appropriate pair of rows (+R0L 
and -RQL, +PCH and -PCH, or +YAW and -YAW) in the response matrix (Table 10, 

11, 12, 13, or 14) are multiplied by the acceleration factors 

A = a/a (24) 

a — 

and 

A e = 6/6 , (25) 

where s represents the steady-state angular acceleration magnitude produced 
by the reinforcing jets. The factor A^ is applied to the row which represents 
a commanded rotation of the same sign as the natural acceleration component, 
and A^ is applied to the row which represents a commanded rotation of the 
opposi te sign. 

To calculate the corresponding propellant consumption rates, columns 
7-12 are multiplied by the factors 

P = A (T/T*) (26) 

a a 
and 

: p e = A a C t/t *) ■ < 27 > 

where, as indicated by Equation (17), T/T* represents the ratio of the elec- 
trical width to the effective width of the average control pulse. Equations 
(26) and (27) represent an implicit variation of specific impulse (l sp ) with 
electrical pulse width. This variation is illustrated for the PRCS jets in 
Figure 4, where the implicitly-defined HFRMP values are superimposed on a curve 
reproduced from Reference 7. 
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EFFECTIVE NOMINAL I CD (SEC) 



Figure 4. Primary RCS Specific Impulse as a Function of Electrical Pulse Width 
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The HFRMP performs the calculations defined by Equations (19) through 
(27), for each one of the Orbiter's three body axes, at every point in the 
trajectory calculation process where state variable derivatives are computed. 
The attitude-control propellant consumption rates and the associated components 
of linear and angular acceleration are summed with their counterparts which are 
produced by other effects (gravity, aerodynamics, and commanded transl ati onal 
thrust), and then integrated by means of a 4th order Runge-Kutta scheme to 
calculate the Orbiter's state vector and propellant expenditure as a function 

of time. 

As a final note of interest with respect to the simplified rate-control 
logic, Equations (19) through (23) indicate that the average acceleration mag- 
nitude in each direction of rotation about a given axis is given by 

a = & = a 2 (C-.014) 2 /16D (28) 


whenever u is of negligible magnitude. When this result is substituted into 
Equations (26) and (27), the average propellant consumption rate for attitude 
control about the axis under consideration is found to be 


_ ctC (C- .014) 
1 6D 


(W + 


V 


(29) 


where W^ and W^, respectively, are the steady-state propellant consumption rates 
for the thrusters that oppose and reinforce the natural acceleration. Further- 
more, the cited equations indicate that the average consumption rate decreases 
to a value of 



(C-.014) 

160 


W 

a 


(30) 


when the natural acceleration magnitude increases to the value 


v = [«(C-.014)] 2 /16D 


(31) 
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and that any further increase in the value of v causes W_ to increase 
again. Equations (29), (30), and (31) very succinctly illustrate the basic 
relationships between the DAP cycle time C, the deadband width D, and the 
average propellant consumption rate due to attitude "deadbanding". The impli- 
cations of these equations are of course subject to the limitations of the 
simplified rate-control logic. One of the more significant limitations is the 
assumption that the natural acceleration is constant everywhere within the 
attitude deadband. This clearly will not be true when the deadband is very 
wide. 


5. COMPARISON OF HFRMP AND SSFS DATA 


For the purpose of verifying the HFRMP OMS/RCS/DAP models, a consider- 
able amount of data were generated with the HFRMP to be compared with similar 
data which were generated by the Space Shuttle Functional Simulator (SSFS) 
and published in References 9 and 10. In the SSFS, thruster firing commands 
are computed at the same frequency as that of the real -world DAP, using a 
jet-select logic and an idealized limit cycle control law representative of 
the Shuttle onboard software. As pointed out previously, the HFRMP does not 
simulate individual RCS thruster firing pulses, but instead calculates the 
average accelerations and propellant rates that result from such pulses in 
an attempt to simulate their long-term effects. Aside from this basic dif- 
ference in simulation technique, the following differences are known to exist 
between the environmental and Shuttle system models that were used in the 
two programs: 

1. The primary reaction control system (PRGS) thrust and moment data 
contained in Tables 1-3 differ significantly in some instances from 
those which were being used in the SSFS when the data were gener- 
ated for References 9 and 10. 

2. Although both programs use the 1962 Standard atmosphere model to 
calculate density as a function of altitude, altitudes are calculated 
with respect to a spherical earth in the SSFS. In the HFRMP, alti- 
tudes are referenced to the Fischer ellipsoid. 

3. The atmosphere is assumed to be stationary (in a geocentric inertial 
frame) in the SSFS, while in the HFRMP it is assumed to rotate with 
the earth. 

4. The aerodynamic coefficients used in the SSFS (at least during the 
generation of data for Reference 9) were representati ve of a doors- 
closed Urbiter configuration. In the HFRMP, aerodynamic coefficients 
are calculated from curve-fit equations that include analytical esti- 
mates of the effect of open doors. 

5.1 PROPELLANT CONSUMPTION : f 

5.1 .1 Translational and Rotat ional Maneuvers 

Tables 16 and 17 contain PRCS propellant requirements for different types 
of translational and rotational maneuvers, as calculated by the SSFS and by 
the HFRMP. Two different cross coupling compensation options (ROT and FULL) 
were used in the HFRMP calculations, and the results obtained from both options 
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Table 16. PRCS Propellant (Lb) Required to Produce a Linear 
Velocity Increment of One Foot per Second 



COMMANDED DIRECTION OF TRANSLATION 

+X 

-X 

+Y 

+1 

-Z 

SSFS (REF. 9) 

27 

30 

42 

23 

43 

HFRMP (ROT COMP.) 

El 

H 

Q 


34 

HFRMP (FULL COMP.) 

28 

34 

55 

25 

m 


(See Table 15 for Orbiter Mass Properties) 


Table 17. PRCS Propellant (lb) Required to Start and Stop 
a Rotational Maneuver Executed at a Rate of 
One Degree per Second 



COMMANDED ROTATION 

ROLL 

PITCH 

YAW 

SSFS (REF. 9) 

18 

22 

24 

HFRMP (ROT COMP.) 

© 



HFRMP (FULL COMP.) 

32 

29 

25 


(See Table 15 for Orbiter Mass Properties) 
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are contained in the cited tables. For each maneuver type, the circled HFRMP 
value represents that generated with the compensation option which most nearly 
approximates the type of compensation applied in the SSFS simulations. 

The SSFS data are taken from Reference 9, The attitude deadbands were 
set at 0.5 degrees per axis and the translational cross coupling thresholds 
(which must be exceeded before any compensation is applied) were set at 0.3 
feet per second per axis for the SSFS calculations. The propellant consump- 
tion data were based on simulated linear velocity increments of one foot per 
second, and angular velocity increments of one degree per second. For maneu- 
vers of this magnitude, the translational cross coupling threshold was ex- 
ceeded only in the case of -Z translation. Therefore, the HFRMP compensation 
mode that most closely approximates the SSFS simulation is FULL in the case of 
-Z translation, while in all other cases it is ROT. No threshold values are 
involved in the HFRMP compensation logic; the effects of cross-coupled acceler 
ations are either removed completely and instantaneously, or not at all. 

Considering the differences in simulation methods and in the basic PRCS 
force and moment data, the HFRMP propellant consumption values shown in Tables 
16 and 17 are reasonably close to those generated by the SSFS. OMS propellant 
consumption data generated by the two programs were identical : 19.9 pounds 
for a velocity increment of one foot per sev md, with an Orbi ter- gross weight 
of 200 017 pounds. 

5.1.2 VRCS Attitude Maintenance 

Table 18 contains vernier reaction control system (VRCS) propellant 
consumption rates for five different attitude-hold options and six orbit 
altitudes. These data, taken from Reference 9, were generated by the SSFS 
with a 40 millisecond DAP cycle time and attitude deadbands of 0.1 degrees 
per axis. Table 19 contains comparable data which were generated by the HFRMP 
using the ROT compensation option and the rate-hold control logic described in 
Section 4.3.2. The orbit inclination was set at 28.5 degrees for all HFRMP 
attitude-maintenance runs. 
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Table 18. SSFS Propellant Consumption Rates (Lb/Hr) ^ For VRCS 
Attitude Maintenance, With 40 mSec DAP Cycle 

(Deadband = 0.1 Deg/Axis) 



ORBIT ALTITUDE (NAUTICAL MILES) 

90 

100 

120 

150 

200 

500 

ATTITUDE 

• i 

X POP INERTIAL HOLD 

9 

6 

3 

2 

1 

1 

Y POP INERTIAL HOLD 

3 

3 

3 

3 

2 

2 

Z POP INERTIAL HOLD 

15 

10 

6 

4 

4 

4 

! X POP ZLV (HEADS DOWN) 

23 

13 

6 

2 

1 

1 

Y TOP ZLV (HEADS DOWN) 

1 

1 

1 

1 

1 

1 


(See Table 15 for Orbiter Mass Properties) 


t - Averaged over a one-orbit period 






























Table 19 


. HFRMP Propellant Consumption Rates (Lb/Hr)' For V RCS 
Attitude Maintenance, With 40 mSec DAP Cycle 

(Deadband = 0.1 Deg/Axis) 



ORBIT ALTITUDE (NAUTICAL MILES) 

90 

100 

120 

150 

mm 

500 

ATTITUDE 

X POP INERTIAL HOLD 

14 

7 

3 

2 

2 

T 

Y POP INERTIAL HOLD 

3 

3 

3 

3 

2 

2 

Z POP INERTIAL HOLD 

17 

11 

8 

6 

6 

4 

X POP ZLV (HEADS DOWN) 

31 

17 

6 

2 

1 

.4 

Y POP ZLV (HEADS DOWN) 

.7 

.5 

.4 

.4 

.4 

.3 


(See Table 15 for Orbiter Mass Properties) 


t - Averaged over a one-orbit period 




































The agreement between the VRCS propellant consumption data produced by 
the two programs is considered to be excellent except at the two lowest altitudes, 
where the differences in the atmosphere and aerodynamic models are sufficient 
to produce significant differences in propellant requirements for the X POP ZLV 
(Orbiter _X axis Perpendicular to the Orbit Plane, Z_axis continuously aligned 
with Local Vertical), the Z POP Inertial Hold, and the X POP Inertial Hold modes. 

The data contained in Tables 18 and 19 represent hourly rates, averaged 
over one orbit period. The presentation of the data in this manner should not 
be taken to mean that propellant consumption is a linear function of time in 
all cases. The HFRMP data shown in Figure 5 illustrates that the instantaneous 
rate can vary significantly over the time span of an orbit period. 

5.1.3 PROS Attitude Maintenance 

Tables 20 and 21 are analogous to 18 and 19, the differences being that 
attitude control was executed by the PROS rather than the VRCS, and the dead- 
bands were widened to 0.5 degrees per axis. The agreement between the results 
of the two programs is not as good as that previously demonstrated for VRCS 
control. The poorer agreement can be attributed partially to the differences 
in the basic PRCS force and moment tables. However, these differences are not 
believed great enough to explain all of the disparities between Tables 20 and 
21, which probably for the most part reflect the inherent limitation of the 
simplified HFRMP logic when it comes to modeling the effects of the relatively 
high degree of overcontrol that often attends the use of the PRCS for attitude 
maintenance. Nevertheless, the agreement is considered adequate to satisfy 
the purposes of the HFRMP' 1 ', especially when taking into account the limitations 
of the computing equipment for which it was designed. 


+i.e., relative motion studies relating to Orbiter/IUS proximity operations 
of comparatively short duration. 
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Table 20. SSFS Propellant Consumption Rates (Lb/Hr) + For PRCS 
Attitude Maintenance, With 40 mSec DAP Cycle 

(Deadband = 0.5 Deg/Axis) 



ORBIT ALTITUDE (NAUTICAL MILES) j 

90 

100 

120 

150 

200 

500 

ATTITUDE 

X POP INERTIAL HOLD 

32 

30 

26 

23 

21 

20 

Y POP INERTIAL HOLD 

22 

16 

12 

10 

10 

18 

Z POP INERTIAL HOLD 

36 

28 

26 

24 

23 

24 

X POP ZLV (HEADS DOWN) 

49 

38 

27 

23 

20 

14 

Y POP ZLV (HEADS DOWN) 

18 

18 

18 

18 

18 

17 


(See Table 15 for Orbiter Mass Properties) 
t - Averaged over a one-orbit period 



Table 21. HFRMP Propellant Consumption Rates (Lb/Hr) f For PRCS 
Attitude Maintenances With 40 mSec DAP Cycle 

(Deadband = 0.5 Deg/Axis) 



ORBIT ALTITUDE (NAUTICAL MILES) 

90 

100 

120 

150 

200 

500 

ATTITUDE 

X POP INERTIAL HOLD 

29 

26 

24 

24 

• 25 

26 

V POP INERTIAL HOLD 

20 

19 

20 

20 

20 

20 

Z POP INERTIAL HOLD 

28 

23 

18 

17 

17 

16 

X POP ZLV (HEADS DOWN) 

42 

34 

25 

25 

26 

27 

Y PDF ZLV (HEADS DOWN) 

27 

27 

27 

. 

27 

27 

. 

27 


(See Table 15 for Orbiter Mass Properties) 

t - Averaged over a one-orbit period 



5.2 TRAJECTORY PERTURBATIONS DUE TO RCS TRANSLATIONAL CROSS COUPLING 


Tables 22 and 23 contain HFRMP data relating to downtrack displacement 
of the Orbiter which results from uncompensated translational cross coupling 
during extended intervals of DAP attitude maintenance. The data were obtained 
by differencing the results of HFRMP runs in which ROT compensation was speci- 
fied (in conjunction with both the V and the P jet select options) with similar 
runs in which the FULL compensation mode was specified in conjunction with the 
P jet select option. The displacements obtained with ROT compensation in- 
cluded the effects of aerodynamic drag and translational cross coupling, while 
the displacements obtained with FULL compensation included the effects of drag 
alone. Since the attitude histories were identical for all three combinations 
of jet-select option and compensation mode, it was reasoned that the displace- 
ments due to drag were essentially the same. It follows that when the displace 
merits obtained with FULL compensation are subtracted from their counterparts 
obtained with ROT compensation, the result represents the displacement due to 
translational cross coupling alone. 

Reference 9 contains no data of the type under consideration here. How- 
ever, Reference 10 contains comparable data for two cases: Y POP ZLV at 150 
miles altitude, with VRCS and PROS control. The Reference 10 data, which 
were derived from SSFS simulations, are enclosed by parentheses in Tables 22 
and 23. The comparisons, although severely limited in number, tend to verify 
the HFRMP results. 

Trajectory deviations due to RCS cross coupling are of considerable 
import in at least two areas of on-orbit activity: 

1. during stationkeeping, where they tend to drive the Orbiter away 
from its desired position relative to a payload, and 

2. during orbit determination, where they tend to pollute navigational 
measurements. 

By no means should the data in Tables 22 and 23 be taken as comprehensive 
general representations of cross coupling effects on the Orbiter ' s trajectory.) 
In the first place, significant displacements can occur in directions other : 
than downtrack. Secondly, the displacements in general are very nonlinear 
functions of time; therefore, the hourly “rates" shown in the cited tables 
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ATTITUDE 


Table 22. HFRMP Downtrack Displacement Rates (Ft/Hr) 1 Due To VRCS 
Translational Cross Coupling, With 40 mSec DAP Cycle 

(Deadband = 0.1 Deg/ Axis) 


X POP INERTIAL HOLD 


Y POP INERTIAL HOLD 


Z POP INERTIAL HOLD 


X POP ZLV (HEADS DOWN) 


Y POP ZLV (HEADS DOWN) 


90 

! 

100 

120 

1 622 

804 

253 

343 

404 

443 

-149 

-82 

-31 

2280 

1214 

447 

34 

27 

23 


ORBIT ALTITUDE (NAUTICAL MILES) 


150 


(See Table 15 for Orbiter Mass Properties) 


200 

500 

254 

234 

469 

423 

0 

V 

14 

45 

19 

21 

19 


Averaged over a one-orbit period 
Data from Reference 10, based on SSFS Simulation 






























Table 23. HFRMP Downtrack Displacement Rates (Ft/Hr)' Due To PRCS 
Translational Cross Coupling, With 40 mSec DAP Cycle 

(Deadband = 0.5 Deg/Axis) 



ORBIT ALTITUDE (NAUTICAL MILES) 

90 

100 

120 

150 

200 

500 

ATTITUDE 

X POP INERTIAL HOLD 

-370 

-249 

-216 

-209 

-206 

-233 

Y POP INERTIAL HOLD 

-150 

-162 

-173 

-190 

-224 

-257 

Z POP INERTIAL HOLD 

-246 

-108 

-24 

-32 

-38 

-41 

X POP ZLV (HEADS DOWN) 

-1137 

-614 

-236 

-no 

-82 

-85 

Y PDF ZLV (HEADS DOWN) 

-823 

-826 

-830 

r (-744)i+ 
-833 

-837 

-937 


(See Table 15 for Orbiter Mass Properties) 

t - Avaeraged over a one-orbit period 

-f - Data from Reference 10, based on SSFS simulation 



cannot be used to estimate even the downtrack displacement alone for any time 
interval appreciably different from the appropriate orbit period"*'. Figures 6, 
7, and 8 illustrate the nonlinearity of the trajectory deviations and the 
diverse nature of the displacement histories for three representative attitude 
maintenance situations. 


The word "rate" is a misnomer for the data under consideration. What the tab! 
actually represent are the displacements at the end of one orbit period, 
normalized to a common interval of one hour to facilitate comparisons between 
orbits having different periods. 
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Figure 6. Orfcriter Displacement Due To Cross Coupling v/ith PRCS Holding Y POP ZLV (Heads Down) Attitude 
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Figure 8. Orb iter Displacement Due To Cross Coupling v/ith VRCS Holding Y POP Inertial Attitude 




6. EFFECTS OF LENGTHENING THE DAP CYCLE TIME 


The Shuttle DAP cycle time has recently been lengthened from the old 
value of 40 milliseconds to a new value of 80 milliseconds (Reference 7) to 
alleviate problems associated with pressure transients in the RCS propellant 
feed lines. This section addresses the influence of that change on trajectory 
perturbations and attitude-maintenance propellant requirements, as determined 
by the results of HFRMP calculations. 

6.1 ATTITUDE MAINTENANCE PROPELLANT COSTS 

Tables 24 and 25, when compared with Tables 19 and 21, indicate that the 
longer DAP cycle time will have a negligible effect on the VRCS propellant 
requirement for attitude maintenance, but that the propellant requirement for 
PRCS attitude maintenance within 0.5 degree deadbands will increase by as much 
as a factor of 5 in some cases. This result can be explained as follows. Use 
of the PRCS for attitude maintenance generally produces a significant degree 
of overcontrol (see Figure 3). In a case of extreme overcontrol , the pulse 
frequency is directly proportional to the minimum effective pulse width, which 
increases by a factor of (.080 - . 01 4 )/ ( . 040 - .014) - 2.5 when the minimum 
electrical pulse width is lengthened from 40 to the new value of 80 milliseconds. 
The average propellant consumption rate is proportional to the product of pulse 
frequency and the electrical pulse width, the latter of course having been 
increased by a factor of 80/40 = 2. 

The preceding rationale is summarized in Equation (29), which also indicates 
that propellant consumption is inversely proportional to the deadband width. 
Therefore, if the DAP cycle time is held at 80 milliseconds, a five-fold in- 
crease in the deadband width (from 0.5 degrees to 2.5 degrees per axis) should 
be expected to reduce the propellant consumption rates in Table 25 to the approx- 
imate levels shown in Table 21 . Table 26, which represents a four-fold in- 
crease in the deadband width, tends to support this conclusion. The same is 
true of Figure 9, which shows the variation of PRCS propellant consumption for 
deadband widths in the range of 0.5 to 10 degrees, as determined from HFRMP 
simulations of a representative attitude maintenance mode. 



ATTITUDE 


Table 24. HFRMP Propellant Consumption Rates (Lb/Hr) 1 For VRCS 
Attitude Maintenance, With 80 mSec DAP Cycle 


(Deadband = 0.1 Deg/Axis) 



ORBIT ALTITUDE (NAUTICAL MILES) 

90 

100 

120 

150 

200 

500 

X POP INERTIAL HOLD 

14 

7 

3 

2 

2 

1 

Y POP INERTIAL HOLD 

3 

3 

3 

3 

3 

2 

Z POP INERTIAL HOLD 

17 

12 

8 

6 

6 

4 

X POP ZLV (HEADS DOWN) 

31 

17 

6 

2 

1 

.4 

Y TOP ZLV (HEADS DOWN) 

.7 

.6 

•5 

.4 

.4 

.3 


(See Table 15 for Orbiter Mass Properties) 


t - Averaged over a one-orbit period 



Table 25. HFRMP Propellant Consumption Rates (Lb/Hr) For PRCS 
Attitude Maintenance, With 80 mSec DAP Cycle 

(Deadband = 0.5 Deg/Axis) 


■' : ' . ■ ' 

.. : ; : "■ , •:= :: , / . 

— - --- -- — - - - - - 

ORBIT ALTITUDE (NAUTICAL MILES) 

90 

100 

120 

— 

150 

200 

500 

ATTITUDE 

X POP INERTIAL HOLD 

119 

126 

1 32 

133 

137 

139 

Y POP INERTIAL HOLD 

110 

111 

113 

115 

115 

112 

Z POP INERTIAL HOLD 

104 

99 

98 

90 

90 

89 

X POP ZLV (HEADS DOWN) 

126 

125 

134 

138 

139 

139 

Y POP ZLV (HEADS DOWN) 

139 

139 

139 

139 

139 

139 


(See Table 15 for Orbiter Mass Properties) 

t - Averaged over a one-orbit period 



ATTITUDE 


Table 26. HFRMP Propellant 
Attitude Maintena 

(Deadband 


X POP INERTIAL HOLD 


Y POP INERTIAL HOLD 


Z POP INERTIAL HOLD 


X POP ZLV (HEADS DOWN) 


Y TOP ZLV (HEADS DOWN) 


90 


33 


26 


30 


44 


34 


(See Table 15 for ( 
Averaged over a one-orbit period 



Consumption Rates (Lb/Hr) + For PRCS 
.nce. With 30 mSec DAP Cycle 

=2.0 Deg/Axis) 


ORBIT ALTITUDE (NAUTICAL MILES) 


100 

120 

150 

200 

500 

31 

31 

32 

33 

34 

26 

26 

27 

27 

26 

25 

22 

21 

21 

21 

37 

30 

33 

34 

34 

34 

34 

34 

34 

35 


3rbiter Mass Properties) 












original; page is 

os ssm suAUxa 


hfrmp data 

150 N. MI. ORBIT ALTITUDE 
V P0P ZLV (HEADS down) attitude hold 


DEADBAND WIDTH, PER AXIS (DEG) 


n? C n P ^°P e i lant Consumption Rate as a Funct 
?yc?e B ?te Width f ° r an 80 - Mi " Second DAP 




6.2 TRAJECTORY PERTURBATIONS 


Table 27, when compared with Table 22, shows that the lengthened DAP 
cycle time has a negligible effect on the magnitude of trajectory deviations 
produced by translational cross coupling arising from VRCS attitude mainte- 
nance. Table 28, on the other hand, indicates (when compared with Table 23) 
that the corresponding effects arising out of PRCS attitude maintenance within 
0.5 degree deadbands are increased by factors greater than 6. The reason that 
the trajectory perturbation factors are higher than the corresponding propel- 
lant consumption factors can be explained, again, in simple terms. Everything 
else being equal, it is logical to expect that the displacements of the Orbiter 
due to uncompensated cross coupling should be proportional to the total momentum 
content of the expelled RCS propellant, which in turn is proportional to the 
product of pulse frequency and the effective width of the thruster pulses (as 
opposed to the electrical width in the case of propellant consumption). It 
has already been pointed out in Section 6.1 that in a severely overcontrolled 
attitude maintenance situation (see Figure 3), the effective pulse width and 
the pulse frequency are each increased by a factor of approximately 2.5 when 
the DAP cycle time is increased from 40 to 80 milliseconds. Therefore, in 

this kind of situation, the trajectory deviations should be expected to in- 

2 

crease by a factor of approximately (2.5) = 6.25. 

As indicated by Table 29, the effect of the lengthened DAP cycle on 
Orbiter trajectory perturbations can be ameliorated by increasing the attitude 
deadband width, in the same manner (and in the same proportion) as it can be 
in the case of propellant consumption; provided of course that other operational 
considerations will permit the use of the wider deadbands. 

Significant increases in the attitude deadband widths may not be feasible, 
for instance, when the Orbiter is stationkeeping at close range with a free- 
flying payload. In such a case, precise translational control of the Orbiter 
must be maintained by the pilot on the basis of the visually observed relative 
motion of the payload. In the presence of large Orbiter attitude excursions, 
it is difficult for the pilot to determine whether a change in the apparent 
position of the payload is actually due to relative translational motion, or 
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ATTITUDE 


Table 27. HFRMP Downtrack Displacement Rates (Ft/Hr) T Due To VRCS 
Translational Cross Coupling, With 80 mSec DAP Cycle 

(Deadband =0.1 Deg/Axis) 


ORBIT ALTITUDE (NAUTICAL MILES) 


X POP INERTIAL HOLD 
Y POP INERTIAL HOLD 


Z POP INERTIAL HOLD 

-149 

-82 

-31 

X POP ZLV (HEADS DOWN) 

2280 

1214 

447 

Y TOP ZLV (HEADS DOWN) 

34 

28 

24 



(See Table 15 for Orbiter Mass Properties) 
t - Averaged over a one-orbit period 




















Table 28. HFRMP Downtrack Displacement Rates (Ft/Hr) ' Due To PROS 
Translational Cross Coupling, With 30 mSec DAP Cycle 

(Deadband =0.5 Deg/Axis) 



ORBIT ALTITUDE (NAUTICAL MILES) 



90 

100 

120 

150 

200 

500 

ATTITUDE 


-1199 

-1247 


-1328 

-1365 

-1547 

X POP INERTIAL HOLD 

-1302 

Y POP INERTIAL -HOLD 

-364 

-881 

-839 

-960 

-1187 

-1463 

Z POP INERTIAL HOLD 

166 

272 

200 

-95 

-221 

-269 

X POP ZLV (HEADS DOWN) 

-1254 

-308 

-569 

-484 

-471 

-528 

Y POP ZLV (HEADS DOWN) 

-5516 

-5512 

-5517 

-5528 

-5544 

-6173 


(See Table 15 for Orbiter Mass Properties) 
T - Averaged over a one-orbit period 




































Table 29, HFRMP Downtrack Displacement Rates (Ft/Hr) : Due To PRCS 
Translational Cross Coupling, With 80 mSec DAP Cycle 

(Deadband = 2.0 Deg/Axis) 



ORBIT ALTITUDE (NAUTICAL MILES) 

90 

100 

120 

150 

200 

500 

ATTITUDE 

X POP INERTIAL HOLD 

-424 

-330 

-334 

-334 


-331 

Y POP INERTIAL HOLD 

-239 

-245 

-247 

-265 

-320 

-385 

2 POP INERTIAL HOLD 

-209 

-70 

-2 

-39 

-58 

-66 

X POP ZLV (HEADS DOWN) 

-1146 

-628 

-249 

-152 

-126 

-134 

Y HP ZLV (HEADS DOWN) 

-1350 

-1352 

-1356 

-1360 

-1365 

-1525 


(See Table 15 for Orbiter Mass Properties) 
t - Averaged over a one-orbit period 


































due to a. rotation of his reference frame (the Orbiter) 1 . On the other hand, 
if the attitude deadbands are not widened, the lengthened DAP cycle may require 
the pilot to execute manual translational corrections much more frequently to 
compensate for the greater translational deviations which attend the automatic 
maintenance of attitude by the DAP. As illustrated in Figures 10 and 11, this 
is especially true if it is necessary to inhibit the upward-firing (+Z) PROS 
jets to minimize plume impingement effects on the payload. 

It is not the intent of this discussion to attempt a general evaluation 
of the overall influence of the lengthened DAP cycle on Orbiter/payload prox- 
imity operations. Such an evaluation is far beyond the scope of this report, 
not to mention the limited capabilities of the HFRMP computer program. The 
intent is rather to point out and to describe some potential operational prob- 
lems that are suggested by results obtained from the relatively simplified 
HFRMP calculations, and which would seem to merit a more comprehensive investi- 
gation involving a real-time man-in-the-loop simulator. 

In this vein, one final observation will be made; this time with regard 
to the potential influence of the lengthened DAP cycle time on tne dynamic 
effects of Orbiter RCS plume impingement on a nearby payload. Even when special 
precautions are taken to fire only those RCS thrusters whose plume centerlines 
extend away from or at right angles to the payload, the plumes expand suffi- 
ciently to cause sensible perturbations of the payload attitude in many cases. 
Again assuming everything else to be equal, the magnitude of these pertur- 
bations can be expected to increase in direct proportion to the total momentum 
content of the expelled RCS propellant. For a given Orbiter deadband, it 
follows then that the lengthened DAP cycle should, increase payload attitude per- 
turb. tions in the same ratio as that previously calculated for Orbiter trans- 
lational perturbations , i . e . , by a factor of approximately 6.25. 


'In some cases, such as in a grappling operation, this makes little difference 
since any kind of relative motion (whether its origin is rotational or trans- 
lational) is objectionable. 
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HFP.MP DATA 
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Figure 10. Hands-Off Stationkeeping Errors v/i t h 40 Millisecond DAP Cycle 






HFRMP DATA 


500 M. MI. ORBIT ALTITUDE 
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ROT COMPENSATION, PZI JET SELECT (UPWARD- 
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V PAYLOAD CG POSITIONS PLOTTED AT ONE-MINUTE 
INTERVALS FOR 5 MINUTES 

/ s' 



Errors with 30 Millisecond DAP Cycle 







7. CONCLUSIONS AND RECOMMENDATIONS 

Comparisons of data generated by the HFRMP and the SSFS computer program 
(Section 5) fail to indicate any gross errors in the HFRMP models of the 
Orbiter OMS/RCS/DAP system. On the basis of these comparisons, it is believed 
that the HFRMP models are sufficiently accurate to produce valid flight pro- 
file designs for proximity operations flight phases. 

The absolute accuracy of the simplified HFRMP attitude control logic 
could not be established for PRCS maneuvers and attitude maintenance because 
the published SSFS data which are currently available were generated prior to 
certain revisions of the PRCS thruster force and moment tables that were in- 
corporated in the HFRMP data base. Recent inquiries directed to personnel 
familiar with the SSFS indicate that the PRCS forces and moments currently used 
by that program are still somewhat different from those used in the HFRMP, by 
reason of differences in the plume impingement effects which are presented in 
Table 3. These inquiries revealed also that a revision of t lie plume impinge- 
ment data is expected in the near future. With this in mind, it was decided 
to postpone any effort toward a more precise correlation of HFRMP and SSFS data 
However, after the plume impingement updates become availab 1 :' and have been 
incorporated in both programs, further data comparisons are recommended for the 
purpose of establishing the accuracy of the HFRMP models more firmly. 

It is believed that the simplified RCS/DAP models utilized in the HFRMP 
are potentially useful in connection with ground-based orbit determination, 
where they might be used to estimate trajectory deviations attributable to 
RCS cross coupling (see Figures 6, 7, and 8). 

Finally, the potential close-range stationkeeping problems that are 
suggested by the HFRMP simulations of 80 mil Ti second DAP cycle effects (Section 
6) are believed to be of sufficient import to merit further investigation by 
means of a real-time man-i n-the-loop simulator; preferably one which models 
the dynamic effects of RCS jet plume impingement on the payload. 
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